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• We present new IRAM Plateau de Bure interferometer (PdBI) 1.3 mm continuum observations at ~1.5" resolution of 

'-O 28 submillimeter galaxies (SMGs), previously discovered with the 870 /im bolometer LABOCA at the APEX tcloscopo 

O from the central 0.7 deg^ of the COSMOS field. Nineteen out of the 28 LABOCA sources were detected with PdBI at 

a > 3cr level of w 1.4 mjy/beam. A combined analysis of this new sample with existing interferometrically identified 
SMGs in the COSMOS field yields the following results: i) > 15%, and possibly up to ~ 40% of single-dish detected 
SMGs consist of multiple sources, ii) statistical analysis of multi-wavelength counterparts to single-dish SMGs shows 
. that only 50% have real radio or IR counterparts, iii) ~ 18% of interferometric SMGs have cither no multi- wavelength 

counterpart or only a radio-counterpart, and iv) ~ 50 — 70% of 2 > 3 SMGs have no radio counterparts (down to an 
1 7-12 /ijy at 1.4 GHz). Using the exact intcrfcrorrictric positions to identify the multi-wavelength counterparts 

allows us to determine accurate photometric redshifts for these sources. The redshift distributions of the combined and 
the individual 1.1 mm and 870 /im selected samples shows a higher mean and a broader width than those derived in 
previous studies. This study finds that on average brighter and/or mm- selected SMGs are located at higher redshifts, 
consistent with previous studies. The mean redshift for the 1.1 mm selected sample {z = 3.1 ± 0.4) is tentatively higher 
than that for the 870 /im selected sample {z = 2.6 ± 0.4). Based on our nearly complete sample of AzTEC 1.1 mm 
SMGs in a 0.15 dcg'^ area, we infer a higher surface density of 2: > 4 SMGs than predicted by current cosmological 
models. In summary, our findings imply that interferometric identifications at (sub-)millimeter wavelengths are crucial 
to build statistically complete and unbiased samples of SMGs. 
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1. Introduction 

1.1. Submillimeter galaxies 

Submillimeter galaxies (SMGs; Ssso^m ^ 5 mJy) are ultra- 
luminous, dusty, starburst galaxies with extreme star for- 
mation rates of order 10'^ Mq yr"-^ (e.g. Blain et al. 2002). 
They trace a phase of the most intense stellar mass build- 
up in galaxies and contribute significantly to the volume- 
averaged cosmic star formation rate density at z = 2 — 3 

20%; Michalowski et al. 2010). Evidence is emerging 
that SMGs represent the progenitors of massive elliptical 
galaxies (e.g. Cimatti et al. 2008; van Dokkum et al. 2008; 
Michalowski et al. 2010), and their enhanced star forma- 
tion properties may be intimately related to the evolution 
of quasi stellar objects (Sanders et al. 1996; Hopkins et al. 
2006; Hayward et al. 2011). 

Spectroscopic and photometric studies of SMGs locate 
them predominantly at redshifts 2-3 (e.g. Chapman et al. 
2005; Wardlow et al. 2011), and only a few z > 4 SMGs 
have recently been detected. Identifying the highest red- 
shift SMGs requires time-consuming, systematic follow-up 
observations to properly identify them against strong lower- 
redshift selection biases. To date ~ 10 z > 4 SMGs have 
been confirmed (Daddi et al. 2009a, 2009b; Gapak et al. 
2008, 2011; Schinnerer et al. 2008; Coppin et al. 2009, 
Knudsen et al. 2010; CariUi et al. 2010, 2011; Riechers et 
al. 2010; Smolcic et al. 2011; Gox et al. 2011; Gombes et al. 
2012). 

Although the number of high redshift (z > 4) SMGs re- 
mains small it appears that their abundance is so high that 
it is only marginally consistent with current galaxy forma- 
tion models (Baugh et al. 2005; Goppin et al. 2009; Smolcic 
et al. 2011). It has also been suggested that high-redshift 
SMGs may be qualitatively different from those at interme- 
diate redshift (Wall et al. 2008). However, such conclusions 
are premature given the significant uncertainties in identi- 
fication of SMGs with optical/IR sources, and therefore the 
appropriate measurement of their redshift. 



1.2. Identifying multi-wavelength counterparts to SMGs 

SMGs are typically first detected with single-dish mm or 
sub-mm telescopes which have a relatively large (10" —35") 
beam size which may include tens of galaxies in the vis- 
ible or NIR. Numerous methods have been applied to 
pinpoint the proper counterparts, such as UV/IR/radio 
star formation indicators or an association with AGN 
(Ivison et al. 2005 2007; Bertoldi et al. 2007; Biggs et 
al. 2010). These techniques are problematic because the 
identification is probabilistic and thus introduces the pos- 
sibility of sample incompletness, contamination, and/or 
bias. Moreover SMGs may be tightly clustered and thus 
blended typical in single dish observations (e.g.. Younger 
et al. 2007, 2009). Interferometric observations of SMGs at 
intermediate-resolution (~ 2") shows that they often do 
not coincide with any galaxy at any wavelength (Younger 
et al. 2009). This may be due to extreme dust extinc- 
tion or due to a very high-redshift for the galaxy. To as- 
sess the overall properties of SMGs including their red- 
shift distribution, it is therefore crucial to follow up single- 
dish detections with high resolution interferometric imag- 
ing. Before the improved sensitivities provided by the up- 
graded IRAM PdBI or ALMA, interferometric follow-up at 



millimeter or submillimeter wavelengths was slow and ex- 
pensive. Only about 50 SMGs have been properly identified 
in various survey fields (Downes et al. 1999, Frayer et al 
2000; Dannerbauer et al. 2002; Downes & Solomon 2003 
Genzel et al. 2003; Kneib et al. 2005; Greve et al. 2005 
Tacconi et al. 2006; Sheth et al. 2004, lono et al. 2006 
Younger et al. 2007, 2009; Aravena et al. 2010a; Ikarashi et 
al. 2011; Tamura et al. 2010; Hatsukade et al. 2010; Wang et 
al. 2011; Chen et al. 2011; Neri et al. 2003; Chapman et al. 
2008; Smolcic et al. 2012). The largest statistically signif- 
icant, signal-to- noise- and flux-limited sample of interfero- 
metrically identified SMGs contains 17 sources drawn from 
the AzTEC/JCMT survey of 0.150° within the COSMOS 
field (Younger et al. 2007, 2009). Here we present PdBI 
observations towards 28 SMGs drawn from the LABOCA- 
COSMOS 0.70° survey (Navarrete et al., in prep.), which 
constitutes the largest interferometric follow-up of SMGs 
drawn from bolometer imaging surveys to date. 

1.3. Determining the redshift of SMGs 

The proper identification with an optical counterpart may 
allow a determination of the SMG redshift through deep 
optical/NIR spectroscopy. Given the ambiguity of identi- 
fications through probability considerations and the opti- 
cal faintness of the counterparts, and the absence of lines 
in particular redshift ranges, this has been a very difficult 
task. The largest SMG sample with spectroscopic redshifts 
to date was established by Chapman et al. (2005), who 
followed-up SMG counterparts identified through deep, in- 
termediate (< 2") resolution radio observations, getting 
redshifts for 76 of 150 targets. 

Where spectroscopic redshifts cannot be measured for 
large samples of SMGs, deep panchromatic surveys such as 
COSMOS or GOODS can measure photometric redshifts, 
which are based on minimization fits of multi-band pho- 
tometry to spectral models (e.g. Ilbert et al. 2009). With 
an optimized choice of spectral models and dense multi- 
wavelength photometric coverage photometric redshifts can 
reach accuracies of a few percent (e.g. Ilbert et al. 2009). 
Although it was not obvious whether common photometric 
redshift templates could be applied to SMG counterparts, 
recent studies confirm that photometric redshifts can be es- 
timated for SMGs, both on statistical and a case by case ba- 
sis (e.g. Daddi et al. 2009a; Wardlow et al. 2010, 2011; Yun 
et al. 2012; Smolcic et al. 2012). Here we further test the 
photometric redshift estimates for SMGs using the largest 
"training set" of SMGs with secure spectroscopic redshifts 
to date from COSMOS. 

In Sec. [2] we describe the data used for our analysis. In 
Sec. [3] we present the PdBI observations towards 28 SMGs 
drawn from the LABOCA-COSMOS survey. In Sec. g] we 
define two samples of SMGs with mm-interferometric de- 
tections in the COSMOS field. Using these in Sec. [5] we 
investigate blending of SMGs, and usually applied statis- 
tical counterpart association methods to single-dish iden- 
tified SMGs. In Sec. 6] we calibrate photometric redshifts 
for SMGs. In Sec. |7| we derive redshift distributions for 
our statistical samples of SMGs with unambiguously deter- 
mined counterparts. We discuss and summarize our results 
in Sec. Hand Sec. ^ We adopt Hq = 70 km s'^ Mpc-\ 
— 0.3, Q\ — 0.7, and use a dust emissivity index of 
/3 — 1, and a Chabrier (2003) initial- mass function if not 
stated otherwise. 



V. Smolcic et al.: mm-imaging of COSMOS SMGs: Redshift distribution 



Table 1. Summary of interferometrically observed COSMOS SMGs besides our work 



Source 


reference 


LABOCA 
source separation Fgroum 
["] [mJy] 


spectroscopic 


redshift 
photometric^ 


mm-to-radio 


1\L _L Jjjvy- 1 


n\ (9^ 


COSLA-60 




12.6 ±3.6 


A t\A 


4-^0-0.20 




A7TFP 9 


("Tl (4) 


COSLA-4 




14.4 ±3.0 


1.125 






A r/TTTT^ Q 
i\Z 1 Hj^y-O 




- 


- 


- 




r on + o.og 

>J.^U_Q 21 






(J-j 


— 


— 


— 








AZ i Il;0-0 




— 


— 


— 


o.y 1 i 


q nt;+0-33 
O.UJ_o 28 




AZ i Il;0-0 










A QOO 
U.OU/ 


n 09+0.13 




AzTEC-7 


(1) 










t-» on + 10 

2.30^^:;^ 




Az i ijjC-8 


(7), (8) 


COSLA-73 




12.3 ±3.6 


3.179 


q 1 7+0.29 
■^■^'-0.22 




AzTEC-9 


(7) 


COSLA-3 




16.4 ±3.3 


1.357 


i-07lg:;j 




AzTEC-10 


(7) 




_ 










AzTEC-11 


(7) 


_ 


_ 


_ 


1.599 


^ rin + 13 




AzTEC-ll-N 


(7) 










1 P-l +0.41 




A ..'IT"/'' 110 


(7) 












> z.oo 


AziEC-12 


(7) 










11/1+0.13 
^■04-0.33 




AzTEC-13 


(7) 


COSLA-158 




11.8 ± 3.9 






> 3.59 


AzTEC-14-E 


(7) 




_ 








> 3.03 


AzTEC-14-W 


(7) 










on + 12 
1-301^:36 




AzTEC-15 


(7) 










3 01 




AzlEO-io 


(9) 








1.505 


i.uy_Q Qg 




JlOOO+0234 


(10) 








4.542 


^■^3_o.o8 




AzTEC/Cl 


(11) 


COSLA-89 




12.4 ±3.7 




5.6 ± 1.2 




Cosbo-1* 


(12) 


COSLA-1 




13.8 ± 1.5 






q Qq+0.68 
•J-O'J_0.49 


Cosbo-3 


(8), (11) 


COSLA-2 




13.1 ± 2.6 


2.490 






Cosbo-8 


(11) 










3.1 ±0.5 




Cosbo-14 


(12) 















Photometric redshifts drawn from the total distribution as described in Sec. [g] and not corrected for any systematic offsets 
*Formany this source is not detected in optical, near- and mid-IR maps/catalogs, therefore we here use the mm-to-radio flux 
based redshift here, which is consistent with the photometric redshift given by Aravena et al. (2010) 

(1) Younger et al. (2007) 

(2) Younger et al. (2009) 

(3) Smolcic et al. (2011) 

(4) Balokovic et al., in prep 

(5) Capak et al. (2010); Riechers et al. (2010) 

(6) Karim et al., in prep 

(7) Younger et al. (2009) 

(8) Riechers et al., in prep. 

(9) Sheth et al., in prep. 

(10) Capak et al. (2009), Schinnerer et al. (2009) 

(11) Smolcic et al., 2012. 

(12) Aravena et al. (2010) 



2. Data 

2.1. The COSMOS Project 

The Cosmic Evolution Survey (COSMOS) is an imag- 
ing and spectroscopic survey of an equatorial 2n° field 
(Scoville et al. 2007). The field has been observed with 
most major space- and ground-based telescopes over most 
of the electromagnetic spectrum. The COSMOS Project 
has obtained very deep broad-band (u*BVgrizJHK) and 
medium and narrow-band imaging data in over 30 optical to 
near-infrared bands. Additionally, there is GALEX, Spitzer 
IRAC/MIPS, Herschel PACS/SPIRE, HST/ACS, XMM- 
Newton, VLA (1.4 GHz and 320 MHz), GMRT (600 & 200 
MHz) data, as well over 25,000 optical spectra (Capak et al. 
2007, Sanders et al. 2007, Scoville et al. 2007, Leauthaud et 



al. 2007, Koekemoer et al. 2009, Frayer et al. 2009, Hasinger 
et al. 2007, Zamojski et al. 2007, Taniguchi et al. 2007; 
Lilly et al. 2007, 2009; Le Floc'h et al. 2009; McCracken et 
al. 2010; Trump et al. 2007, Schinnerer et al. 2004, 2007, 
2010; Smolcic et al., in prep). The inner square degree of 
COSMOS has also been observed in X-rays at a higher res- 
olution and sensitivity with Chandra (Elvis et al. 2009) and 
at mm and submm wavelengths with AzTEC, BOLOCAM, 
MAMBO and LABOCA (Aretxaga et al. 2010; Navarrete 
et al., in prep.; Bertoldi et al. 2007, Aguirre et al., in prep., 
Scott et al. 2008). 

Particularly relevant for the work presented here are the 
deep UltraVista observations of COSMOS which reach 5(t 
(2" aperture AB magnitude) sensitivities of 24.6, 24.7, 23.9, 
and 23.7 in Y, J, H, and Ks bands respectively (McCracken 
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et al. 2012), as well as the VLA 1.4 GHz observations which 
reach a rms of 7-12 /xJy/beam (Schinnerer et al. 2004, 2007, 
2010). We use the updated UV-MIR COSMOS photometric 
catalog (Capak et al. 2007) including all available UV-MIR 
photometric observations. 

2.2. Submillimeter galaxies in the COSMOS field 

The COSMOS field was mapped at mm or submm wave- 
lengths with MAMBO at the IRAM 30m (0.11 deg^; 
1.2mm, 11" angular resolution; Bertoldi et al. 2007), 
BOLOCAM at the CSO (0.27 deg^; 1.1mm; 31" angular 
resolution; Aguirre et al., in prep), AzTEC at the JCMT 
(0.15 deg^; 1.1mm; 18" angular resolution; Scott et al. 
2008), AzTEC at ASTE (0.72 deg^; 1.1 mm; 34" angular 
resolution; Aretxaga et al. 2011), and LABOCA at APEX 
(0.7 deg^; 870 /xm; 27" angular resolution, Navarette et 
al., in prep.). To properly determine the multi- wavelength 
counterparts of the SMGs identified in these surveys, nu- 
merous interferometric and spectroscopic follow-up efforts 
have been made (Younger et al. 2007, 2008, 2009; Capak et 
al. 2008, 2010; Schinnerer et al. 2008; Riechers et al. 2010; 
Aravena et al. 2010a; Smolcic et al. 2011, 2012.; Karim 
et al. 2012, in prep.). To date a sample of 24 interfero- 
metrically identified COSMOS SMGs has been established 
prior to our observations (Table [T]). For 11 of those spec- 
troscopic redshifts are available, either from a dedicated 
COSMOS optical spectroscopic follow-up campaign using 
Keck II/DEIMOS (Capak et al., in prep., Karim et al., in 
prep.), or from CO line observations with mm interferom- 
eters (Schinnerer et al. 2008; Riechers et al. 2010, in prep.; 
Balokovic et al., in prep., Karim et al., in prep., Sheth et 
al., in prep). 



3. PdBI follow-up of LABOCA-COSMOS SMGs 

3.1. Observations 

The COSMOS-LABOCA observation£] reach a rms of 
1.5 mJy per beam (27.6"). The rms increases towards the 
edges of the map, and the catalog was extracted from an 
area of ~ 0.7 deg^ (Navarrete et al. in prep.). Our sample of 
28 COSMOS LABOCA sources selected for PdBI follow-up 
observations (Table [2]) was chosen with a requirement that 
the signal-to-noise m the LABOCA map is S/Ngro^m ^ 
3.8. Eight other LABOCA-COSMOS sources had aheady 
been observed previously with mm-interferometer^ (see 
Table [l]). 

The SMGs in our sample were observed using the PdBI 
during two nights in Oct. /Nov. 2007 (COSLA-10, and 
COSLA-19) and three nights in Oct. /Nov. 2011 (the re- 
maining 26 SMGs) in C- and D-configurations with 6 work- 
ing antennas and the updated PdBI system. All observa- 
tions were done in good/excellent millimeter weather condi- 
tions. During our 2007 observations we used the full 2 GHz 
bandwidth available with the correlator at the PdBI, and 
the receivers were tuned to 232 GHz and 231.5 GHz for ob- 
servations of COSLA-19 and COSLA-10, respectively. The 
total on-source time was 2.3 and 2.2 hrs for COSLA-10, 




^ APEX project IDs: 080.A-3056(A), 082.A-0815(A) and 
086.A-0749(A) 

2 COSLA-1, COSLA-2, COSLA-3, COSLA-4, COSLA-60, 
COSLA-73, COSLA-89, and COSLA-158 



5 10 15 

LABOCA - PdBI separation ["] 

Fig. 1. Distribution of separations between the PdBI 
sources and the corresponding LABOCA-COSMOS 
sources. 



and COSLA-19, respectively. Our 2011 observations were 
done using the WideX correlator covering a bandwidth 
of 3.6 GHz, with receivers tuned to 230 GHz (1.3 mm). 
These observations were performed in snap-shot mode cy- 
cling through the 26 SMGs in each track and observing 
the phase/amplitude calibrator for 2.25 minutes every 19.5 
min. The total on-source time reached is '-^ 43 minutes per 
source. 

Sources J1055-t-018, J1005-I-066, J0923-I-392 were 
used for phase/amplitude calibration, and MWC349, 
J0923+392, 3C84 for flux calibration which we consider ac- 
curate within 10-20%. Calibration and editing was done 
using the GILDAS CLIC package. For each source, the fi- 
nal uv data were collapsed in frequency. The final dirty 
maps reach an rms noise level of 0.55 mJy beam^^ and 
0.39 mJy beam"!, with FWHM beam sizes of 3.3" x 2.3" 
and 3.0" X 2.1" for COSLA-19 and COSLA-10, respectively, 
and an rms of 0.46 mJy beam^^ with FWHM - 1.8" x 1.1" 
for the remaining SMGs. 

3.2. PdBI mm-sources 

We searched for point sources in the dirty 1.3 mm PdBI 
maps within a ~ 14" radius from the phase center, which 
about corresponds to the LABOCA map resolution of 27" . 
Peaks with S/N > 4.5 were considered detections regard- 
less of any multi-wavelength association. When such peaks 
were present in sidelobe-contaminated regions, we tested 
the reality of the sources by cleaning the map by setting a 
CLEAN box around the brightest peak (see Appendix A for 
notes on individual sources). For peaks with 3 < S/N < 4.5 
we required an associated optical, near/mid-IR, or radio 
source within a radius of < 1" . Assuming a Gaussian noise 
distribution, the S/N > 4.5 requirement implies a false de- 
tection rate of ^ 0.15% within a search radius of 14" . A mm 
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Table 2. LABOCA sources observed with the PdBI 
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COSLA-47 


10 


00 


33.1 + 02 26 6.9 


9.0 ± 2.8 


AzTEC/C80 


13.4 


4.1 ± 1.1 








COSLA-48 


10 


00 


24.7 + 02 17 42.3 


6.1 ± 1.7 


AzTEC/C160 


10.0 


3.1 ± 1.2 


Cosbo-7 


12.0 


5.00 ± 0.9 


COSLA-50 


10 


00 


19.0 + 02 16 54.0 


5.6 ± 1.6 














COSLA-51 


10 


00 


11.5 + 02 12 7.1 


6.2 ± 1.7 














COSLA-54 


09 


58 


38.3 + 02 14 2.5 


11.6 ± 4.1 


AzTEC/C13 


8.4 










COSLA-62 


10 


01 


53.2 + 02 20 9.5 


12.5 ± 3.6 














COSLA-128 


10 


01 


38.3 + 02 23 36.1 


11 ± 3.5 














COSLA-161 


10 


00 


15.6 + 02 12 36.0 


5.2 ± 1.7 


AzTEC/C158 


15.6 


3.2+;- 


Cosbo-13S 


5.9 


1.37 ± 0.9 



Table 3. PdBI detections 



source 


position 


Fl.3mm 


S/N 


LABOCA 




redshift 




name 


[J2000] 


[mJy] 




dist. ["] 


S/N870,„„ 


spectroscopic 


photometric^ 


mm-to-radio 


COSLA-5 


10 00 59.521 


-1-02 17 02.57 


2.04 ± 0.49 


4.1 


3.4 


5.0 








COSLA-6-N 


10 01 23.640 


+02 26 08.42 


2.66 ± 0.49 


5.4 


3.4 


4.7 








COSLA-6-S 


10 01 23.570 


+02 26 03.62 


3.08 ± 0.65 


4.8 


7.6 


4.7 




0.48_o 22 




COSLA-8 


10 00 25.550 


+02 15 08.44 


2.65 ± 0.62 


4.2 


6.8 


4.6 




, „o+0.41 
1.83_i 31 




COSLA-9-N 


10 00 13.750 


+01 56 41.54 


1.69 ± 0.47 


3.2 


7.0 


4.5 




„ „„+0.60 
z.oz_2 Q2 




COSLA-9-S 


10 00. 13.829 


+01 56 38.64 


1.87 ± 0.58 


3.2 


5.8 


4.5 




1 90+°^" 




COSLA-ll-N 


10 01 14.260 


-1-01 48 18.86 


2.15 ± 0.62 


3.5 


6.9 


4.4 




n. 7^+0.23 




COSLA-ll-S 


10 01 14.200 


-1-01 48 10.31 


1.43 ± 0.48 


3.0 


2.6 


4.4 




3 00+° " 




COSLA-13 


10 00 31.840 


+02 12 42.81 


2.38 ± 0.61 


3.8 


7.0 


4.3 


2.175 






COSLA-16-N 


10 00 51.585 


+02 33 33.56 


1.39 ± 0.32 


4.3 


3.5 


4.2 




2 16+°-^^ 




COSLA-16-S 


10 00 51.554 


+02 33 32.09 


1.19 ± 0.33 


3.6 


4.3 


4.2 






2 40+0 02 

^■^"-0.51 


COSLA-16-E 


10 00 51.780 


+02 33 33.58 


2.26 ± 0.58 


3.9 


6.0 


4.2 




, „n+3.03 

-i-^^-i.is 




COSLA-17-S 


10 01 36.772 


+02 11 04.87 


3.02 ± 0.57 


5.3 


5.9 


4.2 




70+°-2i 

"■ '"-0.22 




COSLA-17-N 


10 01 36.811 


-1-02 11 09.66 


3.55 ± 0.77 


4.6 


9.1 


4.2 




J-J' _0.22 




COSLA-18 


10 00 43.190 


+02 05 19.17 


2.15 ± 0.48 


4.5 


2.8 


4.2 




2 90+°-^i 

Z.»U_Q 43 




COSLA-19 


10 00 08.226 


+02 11 50.677 


3.17 ± 0.76 


4.1 


11.2 


4.1 






3 98+102 


COSLA-23-N 


10 00 10.161 


+02 13 34.95 


3.42 ± 0.47 


7.3 


1.9 


3.9 




4 00+° '^^ 




COSLA-23-S 


10 00 10.070 


+02 13 26.87 


3.70 ± 0.60 


6.2 


6.4 


3.9 




9 ro + 1-52 
''■00-2.48 




COSLA-33 


10 00 9.580 


+02 19 13.86 


1.78 ± 0.58 


3.1 


6 


3.8 




o 97+0.22 
•'■'^ ' -0.20 




COSLA-35 


10 00 23.651 


+02 21 55.22 


2.15 ± 0.51 


4.2 


3.7 


3.8 




1 91+1™ 

^■■^^-0.64 




COSLA-38 


10 00 12.590 


+02 14 44.31 


8.19 ± 1.85 


4.4 


14.8 


3.7 




2 44+0-12 




COSLA-40 


09 59 25.909 


+02 19 56.40 


3.41 ± 1.02 


3.4 


11.3 


3.7 




, 0(1+0.09 
l^'-^^-O.!! 




COSLA-47 


10 00 33.350 


+02 26 01.66 


3.11 ± 0.59 


5.3 


6.4 


3.6 








COSLA-54 


09 58 37.989 


+02 14 08.52 


3.26 ± 0.65 


5.0 


7.6 


3.6 




2 64+°'^* 




COSLA-128 


10 01 37.990 


+02 23 26.50 


4.50 ± 0.94 


4.8 


10.7 


3.1 




10+01'' 

"■^"-0.00 




COSLA-161 


10 00 16.150 


+02 12 38.27 


2.54 ± 0.74 


3.4 


8.5 


3.1 


0.187 


iq+O-05 
"■^^-0.03 





S/N > 4.5 detections are marked bold-faced 
Photometric redshifts drawn from the total distribution as described in Sec. M and not corrected for any systematic offsets 
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|azTEC 8 



,| AzTEC1 1 | feTE C^ 




AzTEC 



10 

mm flux [mJy] 



Fig. 2. Comparison between LABOCA 870 and PdBI 
1.3 mm fluxes for SMGs (indicated in the panel) detected 
(middle panel) and not detected (top panel) with the PdBI. 
For LABOCA sources identified as multiple PdBI sources 
the individual PdBI source fluxes were added, and for the 
LABOCA sources not detected by PdBI 1.3mm flux es- 
timates for the most likely multi-wavelength counterparts 
were extracted from the PdBI maps (see Sec. |3.3| for de- 
tails). The bottom panel shows the comparison between 
^zTEC 1.1 mm and SMA 890 ^im fluxes (adopted from 
Younger et al. 2007, 2009) for AzTEC/JCMT SMGs in our 
l.lmm-selected sample. The solid line in all panels shows 
the flux ratios for a spectral power law index of 3. 



source association with optical, NIR, MIR, or radio sources 
further decreases the probability that the source is falsfjfj 
Given the surface densities of sources present in various cat- 
alogs the false match probabilities independently estimated 
for each band are 12% (optical), ^ 2% (for each. Ultra Vista 
Y, J, H, Ks, and IRAC 3.6/im), and 0.017% (20 cm radio). 

To further constrain the false match probability, we per- 
formed a source search in the same way as described above, 
but on the inverted, i.e. negative maps. We find only one oc- 
currence of a > 4.5(7 (i.e. A.8a) peak (~ 10" away from the 
phase center and with no multi-wavelength counterpart) 
consistent with the above given false match probability ex- 
pectation. We further find ~ 10% of 3 < S/N < 4.5 peaks 
matched to multi-wavelength counterparts. This suggests 
a ~ 10% false match probability for our 3 < S/N < 4.5 
sources. Hereafter we consider S/N > 4.5 detections as sig- 
nificant and those with 3 < S/N < 4.5 as tentative. 

The 1.3 mm sources and their properties are summa- 
rized in Table [31 Comments on individual sources and their 
PdBI maps and multi-wavelength stamps are presented 
in Appendix A. Of the 28 LABOCA sources observed, 9 
yielded no detection in the 1.3 mm maps (see next Section). 
Six of the 19 detected LABOCA sources break up into mul- 
tiple sources, so that in total we identify 26 submm sources. 
Nine of these 26 sources have S/N > 4.5, 7 have S/N be- 
tween 4 and 4.5, and 10 between 3 and 4. The distribution 
of separations between the LABOCA source position and 
the corresponding PdBI source position is shown in Fig. [l] 
We find a median separation of 6.40" for all sources, and 
5.95" for those with S/Ngyo^m > 3.8. This is consistent with 
the results based on artificial source tests performed on the 
LABOCA map. They result in a positional uncertainty for 
LABOCA sources down to S/NgTo^m = 3.8 of ~ 5.3" with 
an inter-quartile range of 3.1" — 9.8" (Navarrete et al., in 
prep.). 

All detections except COSLA-6-1 and COSLA-6-2 are 
consistent with point-sources at our resolution. We extract 
their fiuxes from the brightest pixel value in the dirty maps. 
The fiux uncertainty is estimated as the rms noise level in 
the map. The fluxes for COSLA-6-1 and COSLA-6-2 were 
obtained by fitting a double Gaussian to the source. All 
fluxes (tabulated in Table 3]) were corrected for the primary 
beam response of the PdBI dishes (assuming a Gaussian 
distribution with HPBW of 21"). 

Scaling the observed 1.3 mm fluxes to the LABOCA 
870 /im, and where available to the AzTEC 1.1 mm, or 
MAMBO 1.2 mm fluxes (Table |2]), yields consistent val- 
ues]^ This is shown in Fig. [2] ana described in more detail 
for each source in Appendix A. This further strengthens 
the validity of our detections. 

3.3. Non-detections 

Nine LABOCA sources remain undetected within our PdBI 
observations. The reasons for this could be that i) the 
LABOCA sources are fainter than our PdBI sensitivity 
limit (Icr ~ 0.46 mJy), ii) the LABOCA sources break 
up into multiple components at 1.5" resolution and are all 



^ If the source is independently detected in various bands then 
the final false match probability is given by the product of the 
individual-band false match probabilities. 

* The fluxes were scaled assuming 5^ oc i/^^'^ where Si, is the 
flux density at frequency v and /3 = 1 the dust emissivity index. 
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fainter than our flux limit or iii) the LABOCA source is 
spurious. To investigate this further we have made use of 
the COSMOS multi-wavelength data by assigning statisti- 
cal counterparts to those LABO CA sou rces given our radio, 
24 /xm, and IRAC data (see Sec. |5.2.1| for details). For 3/9 
sources we find no robust or tentative counterparts while for 
6/9 we find either one or several tentative or robust coun- 
terparts (see Fig. [5] and Appendix A for details) . For the 
latter sources we nave then identified the maximum pixel 
value within a circular annulus of 1" radius in the 1.3 mm 
map. If multiple potential counterparts were present, we 
have summed up the maximum pixel values. Such derived 
1.3 mm fluxes, compared to the LABOCA 870 fj,m fluxes are 
shown in Fig. [2] They agree well with the LABOCA fluxes 
suggesting that the LABOCA sources are not spurious but 
that at interferometric resolution and sensitivity, the source 
is breaking up into multiple-components fainter than our 
1.3 mm sensitivity limit. This is also consistent with the 
results based on artificial source tests performed on the 
LABOCA map which yield that down to a S/Ngyo^^m — 3.8 
5 ± 3 spurious sources are expected (Navarette et al., in 
prep.). 



3.4. Panchromatic properties of PdBI-detected 
LABOCA-COSMOS SMGs 

Twenty-three of the 26 PdBTdetected LABOCA SMGs can 
be associated with multi-wavelength counterparts drawn 
from the deep COSMOS photometric catalog. In addition 
to the UV to MIR photometry from the COSMOS multi- 
wavelength catalog we have added deep YJHK imaging 
from the recent UltraVista Data Release 1. Their pho- 

The COSMOS spec- 



tometry is presented in Table A.l 



troscopic database (Lilly et al. 2007; 2009; Trump et al. 
2007) provided spectroscopic redshifts for the COSLA-13 
and COSLA-161 counterparts. From the 26 SMGs identi- 
fied interferometrically only COSLA-161 was found to be 
associated with X-ray emission in the Chandra-COSMOS 
data (Elvis et al. 2009). 

For each PdBI source we extracted the 1.4 GHz flux 
from the VLA-COSMOS Deep map (Schinnerer et al. 2010) 
using the AIPS task MAXFIT (see Table [A1|). Thirteen of 
the 26 sources 50% with a Poisson error of ±14%) are 
associated with a > 3ct radio peak, where the average rms 
noise level is rmsiAGHz = 9 /iJy/beam. Nine sources are 
detected with S/Ni.4ghz > 4. This radio detection fraction 
does not depend on the significance of the PdBI-source: 
from those with S/N > 4.5 we find five of nine have a radio 
counterpart whereas from those with S/N > 5, three of six 
show a radio counterpart. 



of SMGs with unambiguously determined multi-wavelength 
counterparts. 

In the following we examine two statistically significant 
samples of COSMOS SMGs detected at mm-wavelengths 
at < 2" resolution: 

l.lmm-selected sample: 15 SMGs drawn from the 

1.1 mm AzTEC/JCMT-COSMOS survey at 18" angu- 
lar resolution (AzTEC-1 to AzTEC-15; see Table [T|) 
that form a (S/Ni.imm > 4.5) flux-limited (-Fi.imm^ 

4.2 mJy), 1.1 mm sample. All 15 SMGs were followed- 
up and detected with the SMA at 890 /im, yielding 17 
interferometric sources (as two were found to be mul- 
tiples; Younger et al. 2007, 2009). More details about 
the multi-wavelength photometry of the counterparts 
are provided in Appendix B. 

870^m-selected sample: LABOCA-COSMOS 
sources that were identified at 27" angular resolution 
and confirmed through (sub)mm-interferometry at 
< 2" resolution (Younger et al. 2007; 2009; Aravena 
et al. 2010a; Smolcic et al. 2012., this work). Thirty 
six LABOCA sources were followed-up in total with 
the SMA, CARMA, and PdBI, and 9 resuhed in 
no detection within the PdBI observations down to 
a depth of ~ 0.46 mJy/beam. The remaining 27 
yielded 16 significant (S/N > 4.5) and 18 tentative 
(3 < S/N < 4.5) interferometric (sub)mm-detections. 
For the less significant detections we required an 
association with a source seen at other wavelengths (see 
Table [3] and Table II]). The 16 significant detections 
form the least biasea sample, and we hereafter refer 
to this subsample as the least-biased 870/xm-selected 
sample. 

The sources in the 1.1mm- and 870/im-selected samples 
are summarized in Table [T] and Table [4] Five SMGs belong 
to both samples]^ (see Table |T I . Hereafter we will use these 
two samples to investigate blending, counterpart proper- 
ties, and the redshift distribution of SMGs. For clarity a 
master table of all interferometrically observed SMGs in 
the COSMOS field is given in Table [| 

5. Properties of single-dish detected SMGs when 
mapped at intermediate angular resolution 

In this section we investigate the multiplicity of SMGs de- 
tected at intermediate (< 2") angular resolution, and the 
statistical multi-wavelength counterpart association that is 
commonly applied to single-dish detected SMGs. 



4. Statistical samples of SMGs in the COSMOS 
field identified at intermediate (< 2") resolution 

Our PdBI observations yielded 26 (9 significant S/N > 4.5 
and 17 tentative, 3 < S/N < 4.5) source detections at 
1.3 mm. Combined with previous mm-interferometric de- 
tections of SMGs in the COSMOS field this adds to 50 
SMGs detected with mm-interferometers. To date this is 
the largest interferometric SMC sample. It can be utilized, 
e.g., for a critical assessment of statistical counterpart iden- 
tification methods, and to measure the redshift distribution 



5.1. Blending: Single-dish SMGs breaking-up into multiple 
sources 

In the l.lmm-selected sample of the 15 AzTEC sources 
mapped with the SMA, AzTEC-14 clearly breaks up into 
two sources within the AzTEC beam when observed at ~ 2" 
angular resolution (AzTEC-14-E and AzTEC-14-W), while 
AzTEC- 11 shows extended structure and is best fit by a 
double Gaussian (see Appendix B and Younger et al. 2009 
for details). Thus, in the l.lmm-selected sample only two 

^ AzTEC-l/COSLA-60, AzTEC-2/COSLA-4, AzTEC- 
8/COSLA-73, AzTEC-9/COSLA-3, AzTEC-13/COSLA-158 
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Table 5. Master table of interferometrically observed SMGs in the COSMOS field 



Source 


other names 


(sub)mm-intcrfcroinctry 
observed detected 


statistical 
inter fcrome trie sample 


A'T'TPr' 1 (^'^) 

iT.Z J. JUjO- j. 


COSLA-eO*"', AzTEC/C5*^* SMA, CARMA, PdBI 


V 


1 . Imm , 870/.tiii, least-biased- 8 70/j.rn. 




COSLA-4('*', AzTEC/C3<^' 


SMA, CARMA 


/ 

V 


l.lmm, 870/im, least-biased-870/^m 




AzTEC/C138<^' SMA, CARMA, PdBI 


/ 
V 


1.1mm 


AzTEC-4'- ' 




SMA 


/ 
V 


l.lmm 


AzTEC-S'^'^' 




SMA 


/ 
V 


l.lmm 


AzTEC-6'^'^' 


Az 1 h/O/ClUD^ ^ 


SMA 


, / 
V 


l.lmm 


AzTEC-7'^'^' 




SMA 


. / 

V 


l.lmm 


AzTEC-S'^'^' 


C(JoLA-7o* AziliiC/CJ* ' 


SMA 


,/ 

V 


l.lmm, 870/i.m, least-biased-870//m 


AzTEC-9'^'®' 


COSLA-3^ \ AzTEC/C14^ ' 


SMA 


/ 
V 


l.lmm, 870/im, lcast-biased-870/i.m 


AzTEC-lo'^'^* 




SMA 


,/ 

V 


l.lmm 


AzTEC-11*^'^' 


fRI i-fi! (21 

AzTEC-ll-N'- ^ AzTEC-ll-S^ \ AzTEC/C22'- 


SMA 


, / 
V 


l.lmm 


AzTEC-12*^'^' 


A rx-iTTi 1 1 o ( 2 ) 

AziEO/Clo^ ' 


SMA 


/ 
V 


l.lmm 


AzTEC-13*^'^' 


OUoLA-158^ ' 


SMA 


/ 
V 


l.lmm, 870/im, lcast-biased-870/j.m 


AzTEC-14'^'^' 


AzTEC-14-E*'^\ AzTEC-14-W''^^ 


SMA 


/ 
V 


l.lmm 


AzTEC-15*^'^' 


A 'T~^T7^ 1 1 f 2 ) 

AziEO/ClU^ ' 


SMA 


. / 

V 


l.lmm 


AzTEC-16*^'"'^' 




CARMA 


/ 
V 


~ 


.11000+0234*^'^''' 


A ^'T^TT^ /~l / r'' 1 f7(2) 

AziEC/UiT^ ' 


VLA 


/ 
V 




AzTEC/Cl*^"''' 


/"ly^OT A on(4) 

L-UoLA-^y^ 


A T^ TV T A 
OAKMA 


/ 
V 


870/im 


Cosbo-l*-^' ' 


CUoLA-1^ AziEO/C7^ ' 


SMA 


/ 
V 


870/im, least-biased-870/im 


Cosbo-3^ ' ' 


CUoLA-i^ AziEO/Cb^ ' 


A T> TV ^ A 

CARMA 


V 


870/im, least-biased-870/im 


Cosbo-8* • ' 


AzlEO/Ciio^ 


OAR MA 


V 




Cosbo-14^ ' 




SMA 


/ 
V 




COSLA-5^ ' 


OOSDO-li^^ 


PdBI 


/ 
V 


870/im 


COSLA-6'*'''^' 


COSLA-6-N'^ -', COSLA-6-S^ ^ 


PdBI 


V 


870/im, least-biased-870/im 


COSLA-S'*'''^' 




PdBI 


,/ 

V 


870/im 


COSLA-9' ' ^ 


A „ ("n ■i:r/^//^of2) OTA n m (17) /^i /— . ot a n o(17) 

AziEC/Uo^ OUoEA-y-JN* ', UUbLA-y-b* ' 


PdBI 


/ 
V 


870/im 


COSLA-IO*"^''^' 


L-osbo-D^ 


PdBI 






COSLA-ll*""'''' 


COSLA-11-N'^ '\ COSLA-ll-S*^^^^ 


PdBI 


/ 
V 


870/im 


COSLA-12*''''''' 




PdBI 






COSLA-13**''''' 


AziEL'/L'i4o^ , Oosbo-o^ ' 


PdBI 


,/ 

V 


870/im 


COSLA-14**'''^' 


AzTEC/Cl76^^\ Cosbo-lO*-^^ 


PdBI 






COSLA-16*^''^' 


COSLA-16-N'^ •', COSLA-16-S^ 


PdBI 


,/ 

V 


870/im 


COSLA-17**''''' 


AzlEO/Cli^ , CUbLA-17-iN^ CUdLA-17-o^ ' 


PdBI 


,/ 

V 


870/im, least-biased-870/im 


COSLA-18*^''^' 


AziEo/uyy 


PdBI 


/ 
V 


870/im 


COSLA-19*^''^' 


AzlEO/Oo4* Cosbo-4^ ' 


PdBI 


/ 
V 


870/im 


COSLA-23*''''''' 


/~i ^ 1, „ o(3) i^/^OT A oo /—I o T A no o(17) 

Oosbo-z^ OUoEA-2o-iN^ OUoLA-Jo-o^ ' 


PdBI 


, / 
V 


870/im, least-biased-870/im 


COSLA-25*''''''' 




PdBI 




" 


COSLA-30''*''''' 




PdBI 






COSLA-33*'*''''' 




PdBI 


. / 

V 


870/im 


COSLA-35''*''''' 


AzTEC/C38^^^ 


PdBI 


,/ 

V 


870/im 


COSLA-3S^''''^' 


Cosbo-19^^^ 


PdBI 




870/im 


O(JbLA-40^ ' ' 


AzTEC/C117^^^ 


PdBI 




870/im 


COSLA-47'''''''' 


AzTEC/CSO^^^ 


PdBI 




870/im, least-biased-870/im 


COSLA-48'''''''' 


AzTEC/C160, Cosbo-7'^^ 


PdBI 






COSLA-50*'*''''' 




PdBI 






COSLA-51*'*''''' 




PdBI 






COSLA-54''*''''' 


AzTEC/C13^^^ 


PdBI 


•J 


870/im, least-biased-870/im 


COSLA-62'''''''' 




PdBI 






COSLA-128*''''^' 




PdBI 




870/im, least-biased-870/im 


COSLA-161<''''^' 


AzTEC/Cl58*^\ Cosbo-13S^^^ 


PdBI 




870/im, least-biased-870/im 


(1) Scott ct al. (2008): (2) Aretxaga et al. (2011); (3) Bcrtoldi ct al. (2007); (4) Navarette et al., in prep.: (5) Younger et al. (2007); (6) Younger et al. 
(2009); (7) Smolcic et al. (2011); (8) Balokovic et al., in prep.; (9) Capak ct al. (2010); Riechers et al. (2010); (10) Karim et al., in prep.; (11) Younger 
et al. (2009); (12) Riechers et al., in prep.; (13) Shoth et al., in prep.; (14) Capak et al. (2009); Schinncror et al. (2009); (15) Smolcic et al., 2012; (16) 



Aravena et al. (2010); (17) this work 



of 15 (13% with a Poisson uncertainty of 9%) single-dish 
sources are blended, i.e., they break up into multiple com- 
ponents when observed at intermediate angular resolution. 
The comparison between the single-dish 1.1 mm AzTEC 
and the interferometric SMA 890 //m fluxes for these 15 
sources, shown in the bottom panel of Fig. [2j suggests 
that although the agreement is reasonable, it is possible 
that some faint companions were missed thus potentially 
increasing the fraction of multiples in this sample. 



For the 870^m-selected sample of 36 LABOCA- 
COSMOS SMGs foUowed-up and 27 out of these detected 



with interferometers, 6 SMGfj^ (22% ± 9%) break up into 
multiple sources when observed with interferometers. This 
is within the statistical uncertainties of the results for the 
l.lmm- selected sample. Three more LABOCA SMGs de- 
tected by PdB^ may also consist of multiple components 
(see Appendix A, Aravena et al. 2010b and Smolcic et al. 
2012 for details), and the P-statistics (see next section) 
suggests that at least four of the LABOCA sources not de- 
tected by our PdBI observations!^ are potential blends. This 

® COSLA-6, COSLA-9, COSLA-11, COSLA-16, COSLA-17, 
COSLA-23 
^ COSLA-3, COSLA-5, and COSLA-47 
COSLA-10, COSLA-12, COSLA-48, and COSLA-50 
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Table 4. Statistical samples of SMGs with < 2" angular 
resolution mm-detections in the COSMOS field 



l.lmm-selected 


best 


least-biased- 


best 


sample 


redshift 


87U/im-selectea 

CO T-t-l t-\l p 


reasniit 


AzTEC-1 


4.64+ 


COSLA-1 


'^■"'^ — 0.49 


AzTEC-2 


f.l25+ 


COSLA-2 


2.490+ 


AzTEC-3 


5.299+ 


COSLA-3 


1.357+ 


AzTEC-4 


4.93i;:t?# 


COSLA-4 


1.125+ 


AzTEC-5 


3.97f+ 


COSLA-6-N 


4.oiijir 


AzTEC-6 


0.802+ 


COSLA-6-S 


0.48l«-i«2 


AzTEC-7 


2.30 ±0.10 


COSLA-17-S 


0.701°-^^ 


AzTEC-8 


3.179+ 


COSLA-17-N 


3.54l«-i*# 


AzTEC-9 


1.357+ 


COSLA-23-N 


4.20+°-^^# 
— u.yu 


AzTEC-10 


2 79+J1S 


COSLA-23-S 




AzTEC-fl" 


1.599+ 


COSLA-47 


2.36+°-^* 


AzTEC-lfN** 




COSLA-54 




AzTEC-llS" 


> 2.58 


COSLA-60 


4.64+ 


AzTEC-12 


^■•^^-0.33 


COSLA-73 


3.179+ 


AzTEC-13 


> 3.59* 


COSLA-128 


w-J^W-o.oo 


AzTEC-14-E 


> 3.03* 


COSLA-158 


> 3.59* 


AzTEC-14-W 


i.ou_Q 3e 






AzTEC-15 


n -|y+0.29# 
■^■^ ' -0.37 







Five SMGs belong to both samples; AzTEC-l/COSLA- 
60, AzTEC-2/COSLA-4, AzTEC-8/COSLA-73, AzTEC- 
9/COSLA-3, AzTEC-13/COSLA-158 



** Here we keep the nomenclature given by Younger et al. 
(2009). Note however that AzTEC-ll-S is the northern compo- 
nent of the AzTEC-11 SMG, and AzTEC-ll-N is its southern 
component (see Tab. 1 in Younger et al. 2009) 
+ spectroscopic redshift (see Table 111 for references) 

* mm-to-radio flux ratio based redsmft 

* photometric redshift corrected for the systematic offset of 
0.04(1 4- z), see Fig. [t] with errors drawn from the total x 
distribution 

suggests a fraction of > 6/36 ~ 17%, potentially rising up 
to ~ 40% of LABOCA sources blended within the single- 
dish beam. This is consistent with the fraction obtained if 
only the least-biased-870/im-selected sample is considered 
(see Table |4]). 

5.2. Counterpart assignment methods to single-dish detected 
SMGs 

Here we perform a statistical counterpart assignment for 
the SMGs detected at low angular resolution in our 1.1mm- 
and 870/^m-selected samples, and compare them with the 
exact positions obtained from the interferometers. 

5.2.1. P-statistic 

The most common way to associate single-dish identi- 
fied SMGs with counterparts in higher resolution maps is 
through the P-statistic (Downes et al. 1986), i.e., the cor- 
rected Poisson probability that, e.g., a radio source is iden- 
tified by chance in a background of randomly distributed 
radio/IR sources (Downes et al. 1986; Ivison et al. 2002, 
2005). For a potential radio counterpart of flux density S 
at distance r from the SMG position, P^ = 1 — exp(— Ps[l-|- 
ln(Ps/P3CT)]), where Ps = 1 — exp(— Trr^ns) is the raw prob- 
ability to find a source brighter than S within a distance r 



from the (sub-)mm source, ns is the local density of sources 
brighter than the candidate, and P^a- = Trr^nso- is the criti- 
cal Poisson level, with n^^r being the source surface density 
above the Scr detection level. Robust counterparts are con- 
sidered those with Pc < 0.05, while tentative counterparts 
have 0.05 < Pc < 0.2. 

The commonly used samples search for SMG counter- 
parts are from radio, 24/im, and/or IRAC flux or color- 
selected data (e.g.. Pope et al. 2005, Biggs et al. 2011, Yun 
et al. 2012). The maximum search radius is adjusted to the 
positional uncertainty of the SMG. 

With search radii of 9", and 13.5" for the AzTEC and 
LABOCA SMGs, respectively, we independently computed 
the P-statistics for the potential radio, 24 fim, and IRAC 
color selected {m^ q^„^ — TO4.5^ni ^ 0) counterparts and dis- 
play those in Tables [|] and [t) and Figs.|3)|4) and [6] 

5.2.2. Radio counterparts 

In our l.lmm-selected sample 9/15 (60%) bolometer SMGs 
have radio sources (drawn from the Joint Deep and Large 
radio catalogs with rms ^7—12 /iJy/beam; Schinnerer 
et al. 2007, 2010) within the AzTEC beam. This frac- 
tion is consistent with that found in (sub) mm- surveys (e.g. 
Chapman et al. 2005). Only one mm/SMA source (AzTEC- 
13) in this sample is not associated with a radio source 
present within the single-dish beam (Younger et al. 2009). 
Furthermore, AzTEC-5 and AzTEC-8 each have two P- 
robust radio sources within the AzTEC/JCMT 18" beam. 
In both cases only one radio source is associated with the 
SMA mm-detection. 

Correlating with the Joint VLA-COSMOS Large and 
Deep catalogs, out of the 36 LABOCA SMGs followed- 
up with interferometers, 23 (~ 64%) have radio sources 
(rmsi.4GHz ^7—12 /xJy/beam) within the beam. Of these 
36, 26 were detected at mm-wavelengths with interferome- 
ters (870/im-selected sample), and out of these 26, 17 (65%) 
show radio sources within the LABOCA beam. 

Assigning counterparts to each of these LABOCA 
sources via P-statistic we find that (see Table [6) Fig. [sjind 
Fig. B in 4 cases (2 with S/Ni.3mm > 4.5)0 the ro- 
bust/tentative P-counterpart is not coincident with the 
interferometric-source. Within Poisson uncertainties this is 
consistent with the results from Younger et al. (2009) for 
the l.lmm-selected sample. 

In our 870/im-selected sample COSLA-161 has a mm- 
interferometric detection and two P-robust radio counter- 
parts. Multiple P-tentative radio counterparts are found 
also for COSLA-2, COSLA-5, COSLA-17, and COSLA-73 
(three out of these 5 are significant interferometric detec- 
tions). In all cases, except for COSLA-5, one of the radio 
sources is associated with the inetrferometric-source. 

Combining the above results for our 870/xm-selected 
sample we thus find 4 cases where the robust /tentative 
P-counterpart is not associated with the interferometric 
source, and 4 more ambiguous cases where from the mul- 
tiple robust/tentative P-counterparts found for the SMG 
only one is confirmed by the interferometric source. Taking 
the 26 single-dish SMGs in the 870^m-selected sample this 
amounts to a fraction of 15±8% for the first and latter, sep- 
arately. For the l.lmm-selected sample we find one misiden- 
tified and two ambiguous SMG counterparts assigned via 

^ COSLA-5, COSLA-6, COSLA-8, COSLA-128 
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Fig. 3. 3.6 fim, 24 /im, and 20 cm stamps (30" x30" area) for LABOCA COSMOS sources detected by mm-interferometers 
at < 2" resolution (see Table [3] ) . The bands and sources are indicated in the panels and the names of sources detected 
with interferometers at S/N >4.5 are underlined. The thick yellow circle, 2" in diameter, indicates the mm-interferometer 
position. Robust (square) and tentative (diamond) counterparts determined via P-statistic in each particular (3.6 /im, 
24 /im, and 20 cm) band are also shown (see text for details; see also Table [6]). For each source LABOCA contours in 
Icr steps starting at 2a (with locally determined rms) are overlaid onto the 3T6 /im stamp. 
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Fig. 4. Fig. ^ continued. 



P-statistic. Taking the 15 single-dish SMGs in the 1.1mm- 
selected sample this amounts to 7 ± 7% and 13 ± 7%, re- 
spectively. 



5.2.3. Radio, 24 and IRAC counterparts 



In this section we investigate the agreement between robust 
counterparts determined via P-statistic using radio, 24 /xm, 
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Fig. 5. Same as Fig. but for LABOCA sources not detected within our PdBI observations. 



and IRAC wavelength regimes, and counterparts identified 
via intermediate ^ 2" resolution mm-mapping. 

Where both radio and mid-IR data are available, po- 
tential counterparts to single-dish detected SMGs are com- 
monly selected by searching for P-statistics robust radio 
and 24 /im counterparts. Where no such source can be 
identified, counterparts are searched for among color se- 
lected IRAC sources (ma.gpm — "1,4. > 0). In Table [6] 
and Table [7] we list the resulting P-robust counterparts to 
the LABOCA and AzTEC samples. A summary of the iden- 
tifications is given in Table [8] 

In the 870/im-selected sample we find P-robust coun- 
terparts for 17 out of the 26 PdBI-identified SMGs - irre- 
spective whether these identification are correct or not. In 
total we find 18 P-robust counterparts as COSLA-73 and 
COSLA-161 both have two P-robust counterparts associ- 
ated. From the 18 statistically identified sources, 12 (66%) 
are correct identifications based on our PdBI detections. 
This fraction remains similar if we consider only the single- 
dish detected SMGs with mm-interferometric detections at 
S/N > 4.5, i.e. identified without any prior assumptions 



(i.e. multi- wavelength association) : 11/13 (85%) single-dish 
detected SMGs have P-robust counterparts, in total there 
are 12 P-robust counterparts (as COSLA-73 is in this sub- 
sample) and 7 out of these 12 (58%) match our interfer- 
ometric detections. This amounts to ~ 50% correct iden- 
tifications via P-statistic within the samples analyzed (i.e. 
7/13 for the least-biased- and 12/26 for the 870/xm-selected 
sample). 



In the l.lmm-selected sample we find P-robust coun- 
terparts for 8 of 15 (53%) SMGs with SMA detections 
(Table [tI Fig.|6|. Since AzTEC-5 and AzTEC-8 each have 
two P-robust counterparts, we find 10 P-robust associations 
in total. Seven of the 10 (70%) are coincident with the mm- 
interferometric detections. The fraction remains the same if 
robust anrf tentative statistical counterparts are considered. 
Within the Poisson uncertainties this is consistent with the 
results for the 870/xm-selected sample, i.e. only ~ 50% of 
the single-dish detected SMGs have correct counterparts 
assigned via P-statistic. 
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AzTEC 




Fig. 6. Same as Fig.R but for AzTEC/JCMT/SMA COSMOS sources in our l.lmm-selected sample (see Table [l] and 

Table |7|. The AzTEC/JCMT beam is indicated by the circle in the 3.6 um stamp. 
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Table 6. LABOCA sources observed with nim-interferometers at < 2" resolution, and with counterparts identified via 
P-statistic. Robust statistical P-counterparts [p < 0.05) are marked italic, tentative P-counterparts (0.05 < p < 0.2) are 
shown in regular font, while counterparts identified via mm-interferometry are marked bold-faced. 
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* The radio, MIPS/24/im and IRAC catalogs are available at http://irsa.ipac.caltecli.edu/data/COSMOS/tables/ 



5.3. The biases of assigning counterparts to single-dish 
detected SMGs 

Intensive work has been invested into optimizing techniques 
to determine counterparts to single-dish detected SMGs 
identified at low 10 — 35") angular resolution (e.g. Ivison 
et al. 2002, 2005; Pope et al. 2006; Hainline et al. 2009; 
Yun et al. 2008, 2012). Deep intermediate-resolution ra- 
dio observations, which are less time consuming than sim- 
ilar mm-wave observations, but are expected to trace the 
same physical processes (given the IR-radio correlation; e.g. 



Carilli and Yun 1999; Sargent et al. 2010) have proven ef- 
ficient. However, it was realized that radio-counterpart as- 
signment biases samples to low-redshift (e.g. Chapman et 
al. 2005; Bertoldi et al. 2007). To overcome this, 24 ^m- 
and IRAC color- selected samples have been utilized (e.g. 
Pope et al. 2006; Hainline et al. 2009; Yun et al. 2008). 
Generally such methods identify counterparts to ^ 60% of 
the parent single-dish SMG sample (e.g. Chapman et al. 
2005; Biggs et al. 2011; Yun et al. 2012) and yet the frac- 
tion of misidentifications in these samples remains unclear. 
A further source of bias in such samples is the blending 
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Table 7. AzTEC/JCMT/SMA SMGs with identified robust /tentative counterparts based on the P-statistics. Robust 
statistical P-counterparts {p < 0.05) are marked itahc, tentative P-counterparts (0.05 < p < 0.2) are shown in regular 
font, while counterparts identified via mm-interferometry are marked bold-faced. 
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Table 8. Summary of P-statistic results compared to intermediate resolution mm- mapping 



sample 


radio fraction 
P-statistic* correct ID 


24 /im 

P-statistic* correct ID 


P-statistic* 


m4.5^„ > 
correct ID 


combined 
P-statistic* correct ID 


1 . Imm-sclcctcd'^ 
870/^m-sclcctcd 


8/15 (53.3%) 
11/26 (42.3%) 


7/10 (70%) 
7/12 (58.3%) 


5/15 (33.3%) 5/6 (83.3%) 
6/26 (23.1%) 4/6 (66.7%) 


5/15 (33.3%) 
10/26 (38.5%) 


5/6 (83.3%) 
7/11 (63.6%) 


8/15 (53.3%) 7/10 (70%) 
17/26 (65.4%) 12/18 (66.7%) 



*Only robust (Pc < 0.05) counterparts are considered here 
"Out of 15 bolometer SMGs, 2 robust statistical counterparts are found for each of sources AzTEC-5 and AzTEC-8 (see Table [7| 
''Out of 26 bolometer SMGs, 2 robust statistical counterparts are found for each of sources COSLA-73 and COSLA-161 (see 

Table |6| 



of SMGs within the large single-dish beams. This may po- 
tentially be a severe problem as SMGs have been shown 
to cluster strongly (Blain et al. 2004; Daddi et al. 2009a, b; 
Capak et al. 2010; Aravena et al. 2010b; Hickox et al. 2011) 
and reside in close-pairs (as also suggested by simulations; 
Hayward et al. 2011). Here we provide detailed insight into 
these issues based on statistical samples. 

We generate two unique (870/im- and l.lmm-selected) 
SMG samples with counterparts to LABOCA/APEX and 
AzTEC/JCMT COSMOS SMGs identified via interme- 
diate (< 2" ) resolution mm-mapping (see Sec. |4] and 
Table m. Consistent with results from the literature we 
have found statistical counterparts for ~ 50 — 70% of the 
sources in these samples (see Sec. [5l Table [6] and Table [t]). 
Comparing these with the intermediate (^2") resolution 
mm-detections, we find a ~ 70% match. If there were 
no caveats with the intermediate-resolution mm-detections 
this would imply that statistical counterpart assignment 
methods utilizing deep radio, 24 jiva. and IRAC data (such 
as the one applied here) identify correctly counterparts to 

50% of the parent single-dish samples. Furthermore, it 
would imply that > 15%, and possibly up to ~ 40% of 
single-dish detected SMGs separate into multiple compo- 



nents, with a median separation of ^ 5" , when observed 
at < 2" angular resolution. The misclassification of statis- 
tical assignment is likely intrinsic to the methods applied 
and also due to the break-up of single-dish SMGs into mul- 
tiple components. If indeed a large fraction of SMGs are 
blended within the single dish beams (on scales < 10"), 
this could affect the slope of the (sub)mm counts inferred 
from single-dish surveys as the bright end would be overes- 
timated, while the faint end would be underestimated (see 
Kovacs et al. 2010 for a more detailed discussion). 

We find that radio assignment, relative to near/mid-IR 
wavelength regimes, is the most efficient tracer of single- 
dish detected SMG counterparts (see Table [s]). Thus, as 
already demonstrated by Lindner et al. (2011), who find 
that a 20 cm rms of ~ 2.7 — 5 /iJy identifies radio coun- 
terparts for ~ 90% of SMGs, future deep radio maps with 
EVLA, ASKAP, MeerKAT and SKA will provide efficient 
tracers of SMG counterparts. 

Our samples of LABOCA/APEX and AzTEC/JCMT 
SMGs identified via intermediate (< 2") resolution mm- 
mapping are not complete, but constitute half of the par- 
ent SMG samples (see Scott et al. 2008, Navarette et al., 
in prep.). They are also subject to their own incomplete- 
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Fig. 8. Comparison of spectroscopic and photometric red- 
shifts for 16 starbursts in our COSMOS sample. The pho- 
tometric redshifts were determined using the Michalowski 
spectral templates, and the shown errors are ±3(t errors 
drawn from the distributions of the photometric- redshift 
fits (see Fig. [t] and text for details). The median offset and 
standard deviation of the ^z/{l + Zspoc) distribution are 
indicated in the bottom panel. Note that for z > 3 we find 
a slight, 0.04 • (1 + z), systematic underestimate of the pho- 
tometric redshifts. 



6.1. Calibration and computation of ptiotometric redshifts for 
SMGs 

Photometric redshifts are computed by fitting optimized 
spectral template libraries to the spectral energy distribu- 
tion of a given galaxy, leaving redshift as a free parameter. 
The redshift is then determined via a minimization pro- 
cedure. The quality of the photometric redshifts will depend 
on the choice of the spectral library. To obtain optimal re- 
sults for the population of SMGs using Hyper-z, Smolcic et 
al. (2012) tested three sets of spectral model libraries on a 
sample of eight SMGs in the COSMOS field with counter- 
parts determined via mm-interferometry and with available 
spectroscopic redshifts: 



2T: 



6T: 



nesses and false detection rates within heterogenous data 
sets (assembled from SMA, PdBI, and CARMA observa- 
tions). Thus, although our analysis suggests that roughly 
half of single-dish detected SMGs are correctly identified 
via statistical methods, a more robust insight into these 
issues will have to await further follow-up observations of 
complete samples of single-dish detected SMGs with higher 
sensitivities than the ones presented here and with a uni- 
form rms over the full single-dish beam area. One would 
also preferably want to obtain these data at (at least) two 
separate frequencies. 



6. Distances to submillimeter galaxies 

In this Section we calibrate photometric redshifts for 
SMGs based on a sample of 12 SMGs detected via mm- 
interferometry (and 4 additional high-redshift starburst 
galaxies) in the COSMOS field with spectroscopic redshifts 
spanning a broad redshift range of z ~ 0.1-5.3 (seeTable[l] 
and Table [2| . We optimize the photometric redshift com- 
putation, and apply it thereafter to the remainder of our 
SMC sample. 



Only two - burst and constant star formation history 
- templates drawn from the Bruzual & Chariot (2003) 
library (and provided with Hyper-z). 
Six templates provided by the Hyper-z code: burst, 
four exponentially declining star formation histories 



(star formation rate 



oc e 



'*-/^ where t is time, and 



M: 



0.31, 1, 3 and 5 Gyr) and a constant star forma- 
tion history. This selection of SFH/templates is similar 
to the approach used by Ilbert et al. (2009) to compute 
stellar masses with LePhare. 

Spectral templates developed in GRASIL (Silva et al. 
1998; Iglesias-Paramo et al. 2007) and optimized for 
SMGs by Michalowski et al. (2010). 



They find that all three template libraries yield similar 
results, while the M templates result in the tightest dis- 
tributions. Here we repeat their analysis using a larger sam- 
ple containing 12 SMGs in the COSMOS field with coun- 
terparts determined via mm-interferometry and available 
spectroscopic redshifts. We additionally add to this sam- 
ple 4 sources (Vd-17871, AK03, AK05, AK07), selected in 
the same way as AzTEC-1, AzTEC-5, and JlOOO-l-0234, 
i.e. via criteria identifying high-redshift extreme starbursts 
(Lyman Break Galaxies with weak radio emission; Karim et 
al., in prep.). The photometric redshifts are computed us- 
ing the entire available COSMOS photometry (> 30 bands) 
and the Hyper-z code with a Calzetti et al. (2000) extinc- 
tion law, reddening in the range of Ay = — 5, and allowing 
redshift to vary from to 7. 

The results are shown in Fig. [Tj where we present 
the photometric redshift total distributions for the 16 
sources in our training-set. The overall match between the 
most probable photometric redshift (corresponding to the 
minimum value) and the spectroscopic redshift is good. 
We emphasize that the sample used for this analysis is 
rather heterogeneous in respect of redshift range, detec- 
tions in optical bands, blending, and AGN contribution. 
For example, Cosbo-3 is a blended source not detected in 
images at wavelengths shorter than 1 /im (see Smolcic et 
al. 2012., for details). Constraining its photometric redshift 
well (as shown in Fig. [t]) affirms that our deblending tech- 
niques (described in detail in Smolcic et al. 2012), as well 
as photometric redshift computations work well. Vd-17871 
is a weak SMC (with a CO-line detection, and a contin- 
uum brightness at 1.2 mm of ^ 2.5 mJy; Karim et al., in 
prep.) with substantial AGN contribution identified in the 
IR (Karim et al., in prep). Even in this case our photo- 
metric redshift agrees well with the spectroscopic redshift. 
Note also that within our sample with spectroscopic red- 
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Fig. 7. Photometric redshift total distributions for our SMGs with spectroscopic redshifts. We show results based on 
various sets of spectral models (see text for details): 2T (dotted lines), 6T (dashed-lines) , M (full lines). The spectroscopic 
redshifts are indicated by vertical dashed lines. The source names and the number of degrees of freedom (dof) in the 
photometric redshift minimization are indicated in each panel. The gray-shaded areas in some of the panels indicate 
the redshift range ignored for the determination of the best-fit photometric redshift. The photometric redshift provided by 
the best model (M) and its uncertainty was taken as the minimum value and the 99% confidence interval, respectively, 
both indicated in each panel by the thick and thin red lines. 



shifts there are no catastrophic redshift outliers For two 
sources (AzTEC-3 and AK03) there are two equally proba- 
ble redshift peaks (i.e. Xtot minima). In both cases, however, 
one of those is consistent with the spectroscopic redshift. 
In particular, in the case of AzTEC-3 the low redshift peak 
can be disregarded given that the galaxy is not detected at 
1.4 GHz given the depth of the VLA-COSMOS survey 

In conclusion, comparing the redshift probability dis- 
tributions given the 2T, 6T, and M models, we find that 
the Michalowski (M) models yield the most optimal results 
(i.e. the tightest redshift probability distributions). Hence, 
hereafter we will adopt the Michalowski et al. (2010) spec- 
tral templates for the photometric redshift estimate for 
our SMGs. From the redshift probability distribution for 



The photometric redshift of AzTEC-5 shows the largest de- 
viation from its spectroscopic redshift, but it is still within 2a 
of the (^phot — 2spcc:)/(l -I- Zspcc) distribution (see Fig. |8|. 



a given source we take the most probable redshift (cor- 
responding to that with minimum x^) as the photometric 
redshift of the SMG, and derive the 99% confidence interval 
from its total x^ distribution. The comparison between pho- 
tometric and spectroscopic redshifts is quantified in Fig. [8] 
using the M template library. As already visible from Fig.jT] 
the photometric and spectroscopic redshifts are in very 
good agreement. We find a median of -0.02, and a standard 
deviation of 0.09 in the overall (zphot — 2spcc)/(l + -Zspcc) 
distribution. However, from Fig. |8]it is discernible that the 
systematic offset is higher for higher redshifts. Fitting z < 3 
and z > 3 ranges separately we find a median offset of 0.00, 
and —0.04, respectively, and a standard deviation of 0.11, 
and 0.06, respectively. For comparison, a similar median 
systematic offset (-0.023) has been found by Wardlow et 
al. (2011) for their full sample of LESS SMGs with statis- 
tically assigned counterparts. Yun et al. (2012) find a zero 
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offset for GOODS-South SMGs with statistieally assigned 
counterparts, however their results suggest a slight system- 
atic underestimate of z > 3 photometric redshifts (see their 
Fig. 2), consistent with the results presented here. This sug- 
gests that spectral models used for photometric-redshift es- 
timates could be better optimized for the high-redshift end. 
This is however beyond the scope of this paper, and here 
we will correct the {z > 3) photometric redshifts computed 
for our SMGs for this systematic offset. 

Using the same approach as described above we com- 
pute photometric redshifts for aU SMGs in the COSMOS 
field with multi-wavelength counterparts determined via 
mm-interferometry mapping and without spectroscopic 
redshifts. We present their photometric redshift total 
dist ribu tions (prior to any systematic correction) in Figs, m 
and 10 and tabulate their photometric redshifts (not cor- 
rectea for the systematic offset) in Table [l] and Table |3] 



0.1-10 keV) of (7.5 ± 5) x lO^i, (1.9 ± 1.3) x 10^^^ and 
(3.7 ± 2.5) X lO"''^ erg s~^ (given the marginal detection, 
these values should be considered as upper limits). 

The inferred X-ray luminosities are typical for normal 
galaxies rather than strong AGN {Lx > 10^^ erg s~^; 
e.g. Brusa et al. 2007). This rules out a major AGN con- 
tribution within our SMG sample (consistent with previ- 
ous studies of SMGs; Alexander et al. 2005; Menendez- 
Delmestre et al. 2009), and thus also a significant influence 
of AGN on the accuracy of our photometric redshift esti- 
mates. Furthermore, as demonstrated by the source Vd- 
17871 (see previous Section and Karim et al., in prep.) 
buried AGN only obvious in the IR SED do not appear to 
affect the method. This is consistent with the results from 
Wardlow et al. (2011) who have found that the accuracy of 
photometric redshifts is not affected for SMGs showing an 
IR (8 /im) excess likely due to an AGN component. 



6.2. AGN considerations 

As photometric redshifts are typically computed using li- 
braries for the stellar light only, it may be argued that sub- 
stantial AGN contribution to the UV-MIR SED for some 
SMGs may affect our photometric redshift estimate. Note 
however that only bright Type 1 (broad line) AGN need 
special treatment for photometric redshift estimates (see 
Salvato et al. 2010). For low-luminosity (Seyfert, Type 2) 
AGN, with SEDs dominated by the stellar light of a galaxy 
(e.g. Kauffmann et al. 2003) usual photometric redshift 
computations, as the one presented here, are expected to 
yield satisfactory results. 

To address the AGN issue in our SMG sample we have 
utilized the X-ray data from the Chandra-COSMOS survey 
(Elvis et al. 2009), which provide the most direct way to 
identify AGN associated with the SMGs in our 1.1mm- and 
870/xm-selected samples. 

Only COSLA-161, for which we find a good agree- 
ment between its photometric and spectroscopic redshifts, 
is found to be associated with X-ray emission (note how- 
ever that given the X-ray 0.1-10 keV rest-frame luminosity 
of (6.2 ± 2.7) X 10^° ergs s^^ at the source's low spectro- 
scopic redshift, it is not clear whether the source of X-rays is 
star-formation or emission from the nucleus; see Appendix 
for more details). In order to put further constraints on 
the AGN properties of our SMGs, we derive the average 
X-ray flux in the 0.5-2 keV band using all COSMOS SMGs 
with interferometric positions. This is done in such a way 
that for each SMG we extract the X-ray counts from the 
0.5-2 keV band image within a circular aperture of 1.5" in 
radius, and then convert this to an average X-ray flux. 

For this stacking analysis we only used the so called 
best PSF Chandra mosaic (Elvis et al. 2009), that has a 
continuous coverage of the central 0.5 deg^ of COSMOS 
at 50ks depth, in order to be able to use a small extrac- 
tion region and therefore reduce contamination. The back- 
ground counts were estimated using the stowed Chandra 
background data after normalizing the background image 
to the average background rate in a source-free zone. After 
background subtraction we find a marginal detection at 
a 1.5cr level in the stack with i^o.5-2kcV = (0.9 ± 0.6) x 
10"" erg s"i cm-2. For SMGs at redshifts z =2, 3, and 
4, and assuming a power-law X-ray spectrum with photon 
index 1.8 (typical for AGN), the obtained average flux con- 
verts to average bolometric X-ray luminosities (rest-frame 



7. Redshift distribution of SMGs in the COSMOS 
field 

In this section we present the redshift distributions for our 
1.1mm- and 870^m-selected samples. To derive the redshift 
distributions, we take spectroscopic redshifts if available, 
and otherwise photometric redshifts based on Michalowski 
et al. (2010) spectral templates, and corrected for the sys- 
tematic offset as discussed in the previous Section (see 
Table |4]). 

7.1. Redshift distribution of AzTEC/JCMT SMGs with 
mm-interferometric positions 

Our l.lmm-selected sample contains 17 SMGt[^ with ac- 
curate positions from 890 /im interferometric observa- 
tions at intermediate resolution (^ 2") with the SMA 
(Younger et al. 2007; 2009). Spectroscopic redshifts, based 
on optical (DEIMOS) and/or CO (CARMA/PdBI) spec- 
troscopic observations, are available for 7 out of the 17 
AzTEC/JCMT/SMA SMGs (see Tablejl]). For 7 of the re- 
maining sources we use photometric redshifts, derived as 
described in Sec. [6) Three sources (AzTEC-llS, AzTEC- 
13 and AzTEC-14E) cannot be associated with multi- 
wavelength counterparts in our deep COSMOS images. 
Thus, for these we use the mm-to-radio flux ratio based red- 
shifts, often utilized for the derivation of distances to SMGs 
(CariUi & Yun 1999, 2000). Consistent with the faintness 
at optical, IR, and radio wavelengths the mm-to-radio flux 
based redshifts suggest z > 3 (see Table |4]) for all three 
sources when the PdBI 1.3 mm fluxes and an Arp 220 tem- 
plate are used (following Aravena et al. 2010a). The red- 
shifts for AzTEC/JCMT COSMOS SMGs are summarized 
in Table H 

The redshift distribution for the 17 AzTEC/JCMT 
SMGs mapped by SMA is shown in the left panel of Fig. TT] 
Given that for three sources we only have lower redshift 
limits, we compute the mean redshift using the statisti- 
cal package ASURV which relies on survival analysis, and 
takes upper/lower limits properly into account (assuming 
that sources with limiting values follow the same distribu- 
tion as the ones well constrained) . We infer a mean redshift 
of 3.06±0.37 for the AzTEC/JCMT/SMA sources. We note 



when AzTEC-11 is treated as two separate sources; see 
Appendix B for details. 
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Fig. 9. Same as Fig.^ but for our AzTEC/JCMT/SMA COSMOS SMGs without spectroscopic redshifts. 



that treating the source AzTEC-11 as a single source yields 
a mean redshift of 3.00 ± 0.38. 



7.2. Redshift distribution of LABOCA-COSMOS SMGs with 
mm-interferometric positions 

Unlike the l.lmm-selected sample, the 870/zm-selected one 
is not strictly limited in flux or signal-to-noise. This is be- 
cause a fraction of S/NgTOpm > 3.8 LABOCA-COSMOS 
SMGs have not been detected by our PdBI observations 
at 1.3 mm. The least biased sample that can be con- 
structed from these detections is a sample of the LABOCA- 
COSMOS sources detected with interferometers at mm- 
wavelengths, but without prior assumptions, such as e.g. 
multi-wavelength counterpart associations (as described in 
Sec. [Sj we required that our S/Ni.3„im < 4.5 PdBI sources 
had to be confirmed by independent multi-wavelength de- 
tections). This yields 16 sources (9 PdBI detected with 
S/N> 4.5, and 7 detected with SMA or CARMApJ in our 
least-biased 870/xm-selected sample, listed in Table [4j It is 
hard to assess the completeness of this sample. However, 
if our detection rate of the LABOCA sources can be con- 
sidered random i.e. devoid of any redshift-biases, and the 
properties of the non-detected LABOCA sources are sim- 
ilar to those of the detected ones, then one can assume 
that this subsample reflects the distribution of the par- 
ent LABOCA SMG sample within the same flux limits. 
Furthermore, note that the 9 LABOCA/PdBI sources de- 
tected with S/Ni.3nim > 4.5 within our least-biased 870^m- 
selected sample can be regarded as a 1.3 mm flux-limited 
sample (i^i.smm ^2.1 mJy; given that we reached an rms 
of ^ 0.46 mJy in our PdBI maps). Hence, if the red- 
shift distribution of the least-biased 870/im-selected sam- 
ple is consistent with that of the PdBI sub-sample, we can 
assume that this reflects the distribution of SMGs with 
-Fi.amm > 2.1 mJy. 



Out of the total of 8 LABOCA-COSMOS sources detected 
by CARMA or SMA (see Table[l}, only 7 were detected without 
priors (AzTEC/Cl, 2phot = 5.6 ± 1.2 was detected at 3.2a with 
the CARMA interferometer, and verified by a coincident 4.4(t 
radio source; Smolcic et al. 2012). 



Five out of the 16 SMGs in our least-biased 870/im- 
selected samplep^ have spectroscopic redshifts. For another 
eight we use photometric redshifts as derived in Sec. [6] and 
for 3 others we use the mm-to-radio flux based redshifts. 
Their redshift distribution, shown in the middle panel of 
Fig. 11 is similar to that of the PdBI subsample. Using 



the ASURV statistical package we infer a mean redshift of 
2.59 ± 0.36, and 2.29 ± 0.48 for the least-biased 870/im- 
selected sample, and the PdBI S/Ni smm > 4.5 subsample, 
respectively. 

Given the faintness of the counterpart of COSLA-17- 
S its photometric redshift (0.7j;Q22) is rather poorly con- 
strained, and discrepant compared to the mm-to-radio 
based one (> 4). Thus it is possible that this source is 
at high redshift and further mm-observations of COSLA- 
17-N are required to affirm the reality of this source (see 
Appendix A for details). Nonetheless, excluding COSLA-17 
N and S from the sample, we obtain consistent mean red- 
shifts for the least-biased 870/L(m-selected sample (2.65 ± 
0.38) and the PdBI S/N > 4.5 subsample within (2.34 ± 
0.55). 

We show the redshift distribution of the joint l.lmm- 
selected sample and least-biased 870/im-selected sample in 



the right panel of Fig. [llj A mean redshift of 2.80 ± 0.2^ 
is found. A comparison with results from literature is given 
in Fig. 12 and discussed in detail in Sec. [Sj 



8. Discussion 

8.1. The redshift distribution of SMGs 



In Fig. 12 we compare the (normalized) redshift dis- 
tribution of the AzTEC/JCMT/SMA (left panel) and 
LABOCA/interferometric (middle panel) COSMOS sam- 
ples, and their joint distribution (right panel), with redshift 
distributions of SMGs derived for other surveys (Chapman 
et al. 2005; Banerji et al. 2011; Wardlow et al. 2011; Yun et 
al. 2012). The redshift distribution derived by Chapman et 



" COSLA-2, COSLA-3, COSLA-4, COSLA-60, COSLA-73 
A consistent mean redshift {z — 2.76 ± 0.28) is found if 



AzTEC-11 is treated as a single source. 
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Fig. 10. Same as Fig.jT) but for our LABOCA/PdBI COSMOS SMGs without spectroscopic redshifts. 



al. (2005) is based on a sample of 76 SMGs drawn from var- z = 1.2 — 1.7 in the Chapman et al. sample we supplement 

ious SCUBA 850 /im surveys with counterparts identified it with 19 SMGs with DEIMOS spectra drawn from Banerji 

via radio sources present within the SCUBA beam, and ob- et al. (2011; see their Tab. 2). We combine these two data 

served with Keck I to obtain optical spectroscopic redshifts sets normalizing each by the observed area (556 arcmin^ 
(see their Tab. 2). To account for the redshift desert at 
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Fig. 11. Redshift distribution of our l.lmm-selected sample (left panel), 870/ini-selected sample (middle panel), and 
the two samples combined (with sources present in both samples counted only once; right panel). In the middle panel we 
also show the redshift distribution of our S/Ni imm > 4.5 PdBI-detected LABOCA sample (hatched histogram). Mean 
redshift values, and corresponding errors obtained using the statistical package ASURV, as well as the number of sources 
in each sample are indicated in each panel. Mean redshifts for every sample distribution are also indicated by the thick 
vertical lines. 



for the Banerji et al. , and 721 arcmin^ for the Chapman 
et al. samples; Chapman, priv. com.). 

The distribution published by Wardlow et al. (2011) is 
based on 74 SMGs drawn from the LESS survey at 870 fim 
that could be assigned robust counterparts based on the 
P-statistic (using radio, 24 fim and IRAC data; Biggs et 
al. 2011). Wardlow et al. derived photometric redshifts for 
these galaxies (see their Tab. 2) accurate to (Jaz/{i+z) — 
0.037. Using the P-statistic to associate counterparts to 
SMGs (although using a slightly modified method to that 
utilized by Biggs et al. 2011) Yun et al. (2012) identified 44 
(robust and tentative) counterparts to SMGs detected with 
AzTEC at 1.1 mm in the GOODS-S field. For 16 sources in 
this sample a spectroscopic redshift is used, for 21 a pho- 
tometric redshift was inferred by Yun et al., and for 7 only 
a mm-to-radio based redshift could be derived (see their 
Tab. 3). 

From Fig. [12] it is immediately obvious that the red- 
shift distribution of the COSMOS SMGs is much broader 
compared to that derived from previous surveys, in which 
the SMG counterparts were identified statistically within 
the large bolometer beam. In particular, significant high- 
redshift {z > 4) and low-redshift {z < 2) ends are present. 
In the 870^m-selected sample we find fivcp]out of 16 SMGs 
(i.e. ~ 30%) at z < 1.5. While the redshifts of two of these 
are spectroscopically confirmed, the photometric redshifts 
for the other three show possible secondary (higher) red- 
shift solutions, which are more consistent with their mm- 
to-radio flux based redshifts. In our l.lmm-selected sample 
we find fouiF] out of 17 (23.5%) SMGs at z < 1.5. The 
redshifts for three of these are spectroscopically confirmed. 
Thus, in total we find roughly 20-30% of SMGs at z < 1.5 
(see left and middle panels in Fig. 12). Such low redshift 
SMGs, present in the combined Chapman et al. (2005) and 
Banerji et al. (2011) sample but interestingly missed in the 
1.1 mm AzTEC-GOODS-S and 870 //m-LESS samples, are 
expected in models of the evolution of infrared galaxies (see 

COSLA-3, COSLA-4, COSLA-6S, COSLA-17S, COSLA-128 
^® AzTEC-2, AzTEC-6, AzTEC-9, AzTEC-14W 



e.g. Fig. 7 in Bethermin et al. 2011), and they could be ex- 
plained by cold dust temperatures in these sources (e.g. 
Greve et al. 2006; Banerji et al. 2011). A more detailed 
analysis of the physical properties of these SMGs will be 
presented in an upcoming publication. 

We find significantly more SMGs at the high-redshift 
end {z > 4) in both our 1.1mm- and least-biased-870/xm- 
selected samples, compared to the other surveys. As dis- 
cussed in detail by Chapman et al. (2005) and Wardlow 
et al. (2011) this is likely due to the low-redshift bias of 
statistical counterpart assignment methods. Using statis- 
tical means to overcome this bias Wardlow et al. (2011) 
estimate that ~ 30% (and at most ~ 45%) of all 
SMGs in their sample are at redshifts z > 3. Our com- 
bined AzTEC/ JCMT/SMA and LABOCA/interferometric 
COSMOS data yield that - 50% of the COSMOS SMGs 
with interferometrically identified counterparts are at z > 
3. Exploring these two samples separately, we find that 

- 50% (i.e. 9/17) of the AzTEC/ JCMT/SMA SMGs, and 

- 40% (i.e. 6/16) of the LABOCA/interferometric SMGs 
have z ^ 3. 

It is possible that the discrepancies between the z > 3 
SMG fractions in these different samples are due to their 
different average flux densities. Namely, past studies have 
suggested the existence of a correlation between SMG 
brightness and redshift, in such a way that the brightest 
SMGs lie at the highest redshift (e.g. Ivison et al. 2002; 
Pope et al. 2005; Biggs et al. 2011). The LESS survey source 
flux limit is i^870/im > 4.4 mJy (Biggs et al. 2011). Assuming 
a power-law of 3 this translates into a limit of 2.2 mJy at 
1.1 mm, and 1.3 mJy at 1.3 mm. The AzTEC/ JCMT/SMA 
COSMOS source flux limit is Fi.inim > 4.2 mJy, thus about 
2 times higher, and the LABOCA/PdBI fimit is > 2.1 mJy, 
thus a factor of 1.6 higher compared to the LESS sample. 
Indeed, we find higher mean redshifts (z==3.1±0.4for the 
l.lmm-selected sample, and z — 2.6 ± 0.4 for the 870/im- 
selected sample) and thus also a higher fraction of high- 
redshift sources compared to the results based on the LESS 
survey (z = 2.5 ± 0.3). This is consistent with the sugges- 
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tions from past studies that on average brighter SMGs are 
at higher redshifts. On the other hand, it has also been sug- 
gested that mm-selected samples lie on average at higher 
redshifts, compared to sub-mm-selected samples (e.g. Yun 
et al. 2012). Although our results are also consistent with 
this hypothesis, a more conclusive answer, disentangling 
these degeneracies, will have to await for deeper mm- and 
sub- mm selected samples with interferometric counterparts 
and accurately determined redshifts. 

8.2. High redshift SMGs 

In our 1.1mm- , and 870/im-selected samples we find 9 (3 of 
which have radio counterparts) and 8 (4 of which have radio 
counterparts) z > 3 SMGs. We find 5-^SMGs at z > 4 in 
our l.lmm-selected sample, and and 3-'f^SMGs at z > 4 in 
our 870/^m-selected sample. This corresponds to ^ 30—50% 
of the l.lmm-selected sample, and ~ 20% of the 870/zm- 
selected sample. As our 870/im-selected sample, is not com- 
plete, we can infer only a lower limit for the z > 4 SMG 
surface density of > 3/0.7 « 4 deg~^. The l.lmm-selected 
sample is however nearly complete at the given 1.1 mm flux 
limit of Fi ini„i > 4.2 mJy, and drawn from a uniform area 
of 0.15n°. Four SMGs (AzTEC-1, 3, 4, and 5) in the l.lmm- 
selected sample are found to be at z > 4 (three of those have 
spectroscopic redshifts; see Table |4|. JlOOO+0234, with a 

1.1 mm flux of 4.8 ± 1.5 mJy (i.e. ^ 3cr and thus not 
included in our l.lmm-selected sample), is also spcctro- 
scopically confirmed to be at z > 4 (Capak et al. 2008, 
Schinnerer et al. 2008). Furthermore, only lower redshift 
hmits are available for AzTEC-llS, 13 and 14E. Thus, these 
three SMGs may possibly also lie at z > 4. Hence, these 
5-8 z > 4 SMGs with i^i.imm > 4.2 mJy in the 0.150° field 
yield a surface density in the range of ~ 34 ± 14 deg~^ to 
~ 54± 18 deg~^ (Poisson errors). Both values are substan- 
tially higher than what is expected in cosmological models 
(Baugh et al. 2005; Swinbank et al. 2008; see also Coppin et 
al. 2009, 2010), even if the AzTEC/JCMT-COSMOS field 
were affected by cosmic variance of to a factor of 3 over- 
density, as suggested by Austermann et al. (2009). 

Based on the galaxy formation model of Baugh et al. 
(2005; top-heavy IMF; A cold dark matter cosmology) a 
surface density of ~ 7 deg~^ for z > 4 SMGs with 850 iim 
fluxes brighter than 5 mJy is expected (see also Swinbank et 
al. 2008; Coppin et al. 2009). As the AzTEC/JCMT/SMA 

4.2 mJy flux limit at 1.1 mm translates to about a factor of 
two higher flux (i.e. 9.6 mJy) at 850 fiJy, the models would 
predict an even lower surface density at this flux threshold. 

In the GOODS-N field to date four z > 4 SMGs were 
found (Daddi et al. 2009a, b; Carilh et al. 2011). Given 
the 10 X 16.5 arcmin^ area (but with a highly non uniform 
rms in the SCUBA map with average la = 3.4 mJy; Pope 
et al. 2005) this implies a surface density of > 87 deg~^, 
an order of magnitude higher than predicted by the mod- 
els. The GOODS-N z > 4 SMGs are however associated 
with a protocluster at z ~ 4.05 which increases the sur- 
face density value. The COSMOS z > 4 SMGs were se- 
lected from a larger field, and although an overall over- 



" AzTEC-1, AzTECS, AzTEC-4, AzTEC-5, J1000-)-0234, and 
possibly AzTEC-llS, AzTEC-13, and AzTEC-14E for which 
only lower redshift limits are available; see Table m and Tabled 
COSLA-6-1, COSLA-23-N, COSLA-60, and possibty 
COSLA-158 with only a lower-redshift limit; see Table ^ 



density of bright SMGs was found in the AzTEC/JCMT- 
COSMOS field (Austermann et al. 2009), the z > 4 SMGs 
do not seem to be associated with each other (Capak et 
al. 2009, 2010; Schinnerer et al. 2009; Riechers et al. 2010; 
Smolcic et al. 2011). We still find significantly more z > 4 
SMGs than current models predict. If the AzTEC/JCMT 
COSMOS SMGs are representative of the overall SMG pop- 
ulation (-Fi.imm > 4.2 mJy), then our results imply that 
current semi-analytic models underpredict the number of 
high-redshift starbursts. 



9. Summary 

We presented PdBI continuum observations at 1.3 mm with 
~ 1.5" angular resolution and an rms noise level of ~ 
0.46 mJy/beam towards 28 SMGs selected from the (single- 
dish) LABOCA-COSMOS survey of 27" angular resolu- 
tion. Nine SMGs remain undetected, while the remainder 
yields 9 highly significant (S/N > 4.5) and 17 tentative 
(3 < S/N < 4.5 with multi-wavelength source association 
required) detections. Combining these with other single- 
dish identified SMGs detected via intermediate (< 2" ) an- 
gular resolution mm-mapping in the COSMOS field we 
present the largest sample of this kind to-date, contain- 
ing 50 sources. Based on 16 interferometrically confirmed 
SMGs with spectroscopic redshifts, we show that photo- 
metric redshifts derived from optical to MIR photometry 
are as accurate for SMGs as for other galaxy populations. 
We derived photometric redshifts for those SMGs in our 
sample which lack spectroscopic redshifts. 

We distinguish two statistical samples within the total 
sample of 50 COSMOS SMGs detected at < 2" angular 
resolution at mm- wavelengths: i) a l.lmm-selected sample, 
forming a significance- (S/Ni imm > 4.5) and flux- limited 
(Fi.inim > 4.2 mJy) sample containing 17 SMGs with inter- 
ferometric positions drawn from the AzTEC/JCMT 0.150° 
COSMOS survey, and ii) a 870/im-selected sample, contain- 
ing 27 single-dish SMGs drawn from the LABOCA 0.70° 
COSMOS survey and detected with various (CARMA, 
SMA, PdBI) mm-interferometers at intermediate angular 
resolution. 

Within our samples we find that > 15%, and up to 
~ 40% of single-dish identified SMGs tend to separate into 
multiple components when observed at intermediate angu- 
lar resolution. 

The common P-statistics counterpart identification cor- 
rectly associates counterparts to ^ 50% of the parent single- 
dish SMG samples analyzed here. 

We derive the redshift distribution of the SMGs with 
secure counterparts identified via intermediate < 2" reso- 
lution mm-observations, and compare this distribution to 
previous estimates that were based on statistically identified 
counterparts. We find a broader redshift distribution with 
a higher abundance of low- and high-redshift SMGs. The 
mean redshift is higher than in previous estimates. This 
may add evidence to previous claims that brighter and/or 
mm-selected SMGs are located at higher redshifts. 

We derive a surface density of z > 4 SMGs (i^i.imm ^ 
4.2 mJy) of ~ 34 — 54 deg~^, which is significantly higher 
than what has been predicted by current galaxy formation 
models. 



22 




Fig. 12. Normalized redshift distributions for SMGs drawn from various studies in the literature, indicated in the panel. 
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Appendix A: Notes on individual 

LABOCA-COSIS/IOS targets observed with the 
PdBI 

Here we present detailed notes on individual LABOCA- 
COSMOS (i.e. COSLA) SMGs observed with PdBI at 
1.3 mm and ^ 1.5" resolution. Wc use a spectral index 
of 3, i.e. assuming oc where is the flux density 

at frequency v and /3 = 1 the dust emissivity index, to con- 
vert fluxes from/to various (sub-)mm wavelengths (if not 
stated otherwise). 

COSLA-5 

COSLA-5 is a S/N = 4.1 detection located at a = 10 
00 59.521, 5 +02 17 02.57. The PdBI-source is found at 
a separation of 3.4" from the LABOCA source center. The 
1.3 mm flux density of the source is 2.04±0.49 mJy. Scaling 
this to the MAMBO (1.2 mm) and LABOCA (870 /xm) 
wavelengths we find 2.6±0.6 and 6.9±1.7 mJy, respectively. 
This is slightly lower than the extracted (and deboosted) 
MAMBO and LABOCA fluxes (4.78±1, and 12.5±2.6 mJy) 
and may indicate the presence of another mm-source, not 
detected in our PdBI map. Based on the rms reached in 
the PdBI observations, we can put a 2>a upper limit to this 
potential second source of 1.4 mJy. 

The PdBI peak is 1.3" away from a source that is inde- 
pendently detected in the optical (z+ = 22.5), UltraVista, 
and IRAC bands. The photometric redshift of this source 
is well constrained, Zphot = 0.85tQ;oe- Interestingly, at a 



separation of 1.1" towards the SW of the PdBI detection 
we find a faint source present only in the Ks band images, 
but not included in the catalogs (it is present in both the 
WIRCam and UltraVista images). 

The mm-source is not associated with a radio de- 
tection suggesting a mm-to-radio fiux based redshift for 
COSLA-5 of 2;mm/radio ^ 3.8. Here we take the opti- 
cal/IRAC/Ultra Vista source as the counterpart noting that 
given the ~ 4it significance of the PdBI peak, and a separa- 
tion of > 1" from multi-wavelength sources, further follow- 
up is required to confirm this source and its redshift. 

COSLA-6 

Two significant (S/N > 4.5) sources are detected within 
the COSLA-6 PdBI map. 

C0SLA-6-N (S/N=5.4, a =10 01 23.64, 5 =+02 26 
08.42) has a 230 GHz (1.3 mm) flux density of 2.7+0.5 mJy. 
No IRAC/UltraVista source is found nearby. The closest 
source is a faint optical (no Ultra Vista/IRAC/radio) source 
2.0" away. Given the high significance of the mm-detcction 
the mm-positional accuracy is ~ 0.3". Thus, it seems un- 
likely that this optical source is the counterpart of the 
mm-detection. C0SLA-6-N is however coincident with a 
2.1(7 peak in the radio map (Fi.4ghz = 19.5 + 9.4 /uJy). 
Based on this, we infer a mm-to-radio flux based redshift 

of .2^mm— radio — 4.01_q gg. 

C0SLA-6-S (S/N = 4.75, a =10 01 23.57, 5 =+02 
26 03.62) has a 1.3 mm flux density of 3.1 + 0.6 mJy. It 
might be associated with a source detected in the opti- 
cal (separation= 0.5"; = 26.15), but not in near- or 
mid-IR. It is coincident with a 3.3ct peak in the radio 
map (i^i.4GHz = 33.3 + 10.1 /xJy). Based on the multi- 
wavelength photometry, we infer a photometric redshift of 
= 0.481 

0.22 ^'^'^ ttii^ source. A second potential red- 
shift solution (although not as likely as the first one) exists 
at 2 ~ 4, and it is supported by the mm-to-radio flux based 
redshift, ^mm/radio = 3.44l^;^^ . 

The combined 1.3 mm fluxes of C0SLA-6-N and 
C0SLA-6-S, scaled using a spectral index of 3, yield an 
expected flux density of 19.4 + 2.7 mJy at 870 ^m. This 
is in very good agreement with the deboosted LABOCA 
870 /mi flux of 16.0 + 3.4. 

CQSLA-8 

COSLA-8 is detected at S/N = 4.2, and located at 
a = 10 00 25.55, 5 = +02 15 08.44. Its 1.3 mm flux density is 
^^i.3mm = 2.65 ± 0.62 mJy. Using a spectral index of 3 this 
extrapolates to an 870 /zm flux density of 8.9 ±2.1 mJy, 
in good agreement with the LABOCA deboosted flux of 
6.9+1.6 mJy. The source is located in a crowded region. A 
radio, IRAC, UltraVista detection is present at a separation 
of ~ 3", however the closest source to the mm-source (sep- 
aration=1.0") is detected only in the optical (i+ = 27.4). 
A 3.3ct peak is found at the mm-position in the radio map 
(^'i.4GHz = 26.2 + 8.0 /uJy). The most probable photometric 
redshift for this source is ^phot = 1.83l°;3;^, however with a 
rather flat distribution as reflected in the uncertainties. 

COSLA-9 

We identify two 3.2(7 peaks at a = 10 00 13.83, 
5 = +01 56 38.64 (C0SLA-9-S; 5.8" away from 
the LABOCA source center), and a = 10 00 13.75, 
5 = +01 56 41.54 (C0SLA-9-N; 7" away from the LABOCA 
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Fig. A.l. Optical to radio stamps for the LABOCA-COSMOS SMGs detected by PdBI. Names of sources with S/N > 4.5 
in PdBI maps are underhned. 
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source center). C0SLA-9-S can be matched to an opti- 
cal/IRAC/UltraVista source (z+ = 24.8, separation=0.8"), 
while the closest source to C0SLA-9-N is detected only in 
the optical {i+ = 26.1; separation=0.4" ) however it is only 
1.3" away from an optical/UltraVista/IRAC source. 

The primary beam corrected 1.3 mm flux densities of 
C0SLA-9-N and C0SLA-9-S are 1.69 ± 0.47 mJy, and 
1.87±0.58 mJy, respectively. Added together, and scaled to 
the LABOCA (13.2 ±2.1 mJy) and AzTEC (6.3±1.0mJy) 
frequencies yields a good match to the deboosted LABOCA 
(14.4 ± 3.3 mJy) and AzTEC (8.7 ± 1.1 mJy) fluxes. Given 
the low significance of the sources, further follow-up is re- 
quired to confirm their reality. 

COSLA-10 

No significant source is present in the 1.3 mm m ap. 
The statistical counterpart association (see Sec. 5.2.1 and 
Fig. [5] for details) suggests three separate potential (tenta- 
tive) counterparts to this LABOCA source. The sum of the 
extracted 1.3 mm fluxes (taken as maximum flux within 
a circular area of 1" in radius centered at the statistical 
counterpart), corrected for the primary beam response, is 
2.43 mJy. This yields a flux of 8.1 mJy when scaled to 
870 fj,m, in very good agreement with the LABOCA flux of 
7.3±1.7 mJy (see Fig.[2|. This suggests that the LABOCA 
source may be fainter at 1.3 mm than can be detected given 
our PdBI sensitivity and that it possibly breaks up into 
multiple components when observed at 1.5" resolution. 

CQSLA-11 

The brightest peaks in the PdBI map are at 3.5(7 
(COSLA-ll-N) and 3a (COSLA-ll-S). 

COSLA-ll-N is located at a = 10 01 14.260, ^ = 01 48 
18.86, and it can be associated with a faint optical detec- 
tion 0.64" away (i+ = 27.75). Its 1.3 mm flux density is 
-Fi.Smm = 2.15±0.62 mJy, and the source is not detected in 



the radio map. Based on the multi-wavelength photometry 
of the counterpart of COSLA-ll-N, we flnd a photomet- 



ric redshift of 



^phot 



= 0.75: 



-0.23 



The mm-to-radio based 
> 3.6. 



redshift however suggests 2;inm/radi( 

COSLA-ll-S (a = 10 01 14.200, S = +01 48 10.31) 
is only 2.6" away from the center of the LABOCA source, 
and it coincides (separation=0.5") with independent opti- 
cal and UltraVista H-band, and UltraVista J-band detec- 
tions. Although fairly low S/N, the UltraVista detection 
increases the probability that it is a real source. Its 1.3 mm 
flux density is Fi.snun — 1.43 ± 0.48 mJy. Our photomet- 
ric redshift estimate yields two, almost equally probable, 
redshifts at z ~ 0.2 and z ~ 3. Given that the SMG is not 
detected in the radio map, the higher redshift solution, also 
consistent with the mm-to-radio flux ratio based redshift 
(■^mm/radio ^ 3), is morc likely. We thus adopt the higher 
redshift solution for this source yielding Zphoto — S-OO^.q'^'^ 
(where the errors reflect the 99% confidence interval derived 
using only z > 1.5 values). 

The combined flux densities of the two detections, scaled 
to 870 /im, yield a flux density of 12.0 ± 2.6 mJy at this 
wavelength. This is in good agreement with the LABOCA 
870 /im flux density of 19.4 ± 4.5 mJy, and thus further 
affirms the reality of the sources. 

COSLA-12 

No significant source is present in the 1. 3 mm map. The 
statistical counterpart association (see Sec. 5.2.1 and Fig. 
for details) suggests two separate potential (one robust am 
one tentative) counterparts to this LABOCA source. The 
sum of the extracted 1.3 mm fluxes (taken as the maxi- 
mum flux within a circular area of 1" in radius centered 
at the statistical counterpart), corrected for the primary 
beam response, is 2.68 mJy. Scaled to 870 /im this implies 
a flux of 8.9 mJy which is lower than the LABOCA flux of 
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18.3±4.2 mJy (see Fig.[2]). This suggests that the LABOCA 
source may be fainter at 1.3 mm than can be detected given 
our PdBI sensitivity and/or it breaks up into multiple com- 
ponents when observed at 1.5" resolution. 

COSLA-13 

COSLA-13 is detected at S/N = 3.9 [a = 10 00 
31.840, S = +02 12 42.81). Its 1.3 mm flux density is 
-Fi.3mm = 1-37 ± 0.61 mJy This SMG is detected within 
the LABOCA (870 A*ni), AzTEC (1.1 mm), and MAMBO 
(1.2 mm) surveys. Scaling the PdBI 1.3 mm flux to 870 /im, 
1.1 mm, 1.2 mm (using a slope of 4.6 which corresponds 
to the mean slope between the AzTEC/LABOCA and 
AzTEC/MAMBO detected fluxes) we find flux densities 
of 8.8 ± 3.9 mJy (870 fim), 2.9 ± 1.3 mJy (1.1 mm), 
and 2.0 ± 0.9 mJy (1.2 mm). These are consistent with 
the AzTEC/LABOCA fluxes, and slightly lower than the 
MAMBO flux (note that we find consistent results when 
using a slope with a spectral index of 3). 

The PdBI mm-source is coincident with an op- 
tical/UltraVista/IRAC/radio source (separation=0.55") 
with an optical spectrum at redshift Zgpcc — 2.175. 

COSLA-14 

No significant source is present in the 1. 3 mm map. The 
statistical counterpart association (see Sec. 5.2.1 and Fig. [5] 
for details) suggests one robust statistical counterpart to 
this LABOCA source. The extracted 1.3 mm flux (taken as 
maximum flux within a circular area of 1" in radius centered 
at the statistical counterpart), corrected for the primary 
beam response, is 1.64 mJy. The scaled 870 /im flux of 
5.5 mJy is fairly consistent with the LABOCA flux of 9.0 ± 
2.1 mJy when scaled to 870 /im (see Fig. [2]). This suggests 
that the LABOCA source may be fainter at 1.3 mm than 
can be detected given our PdBI sensitivity. 

COSLA-16 

A significant extended source is found ^ 3.5" away from 
the LABOCA source center. It is best fit by a double- 
Gaussian (using the AIPS task jmfit and fixing the width of 
the Gaussians), yielding two sources located at a = 10 00 
51.5854, 6 = +02 33 33.5648 (C0SLA-16-N) and a 10 
00 51.5541, S = -1-02 33 32.0948 (C0SLA-16-S). The 2- 
Gaussian fit yields 1.3 mm flux densities of i^i 3 mm = 1.39± 
0.32 mJy (C0SLA-16-N) and Fi 3 mm = 1-19 ± 0.33 mJy 
(C0SLA-16-S). 

C0SLA-16-N can be associated with an opti- 
cal/UltraVista/IRAC/radio source (separation — 0.79"). 
C0SLA-16-S is not associated with a separate source in the 
multi-wavelength catalogs. The radio emission associated 
with the position of C0SLA-16-N is significant (-Fi.4ghz = 
95.6 ± 10.1 /iJy), while a 3.3a peak, that seems to be 
the extension of the significant radio source, is associated 
with C0SLA-16-S (i^i.4GHz = 33.3 ± 10.1 fiJy). The multi- 
wavelength photometry of C0SLA-16-N implies a photo- 
metric redshift of 2;phot = 2.16l^Q25. '^^^ mm-to-radio based 



S'phot 



= 1.25 



-3.03 
-1.15' 



redshift inferred for C0SLA-16-S is 



a/radi( 



= 2.40 



0.62 
-0.51' 



suggesting it is associated with C0SLA-16-N. 

A third 3.9ct peak (C0SLA-16-E) with 1.3 mm flux 
density of F1.3 mm — 2.26 ± 0.58 mJy is found 6" east 
of the LABOCA source center and it is coincident with 
a faint optical source (separation = 0.41", i+ = 29.20). 
Our photometric redshift computation yields s redshift of 



however (as also reflected in the error) 
the distribution is fairly flat below z < A thus mak- 
ing all redshifts below z ~ 4 almost equally probable. The 
mm-to-radio flux based redsfhit suggests ^^mm/radio ^ 3.7. 

The combined 1.3 mm flux density of the 3 identified 
sources, scaled to 870 fim yields 16.3 ± 2.5 mJy, in very 
good agreement with the deboosted LABOCA flux (14.0 ± 
3.6 mJy). 

COSLA-17 

Two significant S/N > 4.5 detections are found within 
the PdBI map. 

C0SLA-17-S is detected at high significance (S/N — 
5.3; a = 10 01 36.772, 6 = +02 11 04.87). Its 1.3 mm flux 
density is i^i.3 mm = 3.0 + 0.6 mJy, and it can be associated 
with a faint source 0.23" away with mNBsie = 26.2. No IR 
or radio source is associated with this detection. We find 
a photometric redshift for this source of Zphot — 0.7^0^2 , 
while the mm-to-radio based redshift suggests ^mm/radio ^ 
4. 

C0SLA-17-N is found at S/N=4.6 (a = 10 01 36.811, 
6 — +02 11 09.66) with a 1.3 mm flux density of 
-Fi.Smm = 3.55 + 0.677 mJy. It is perfectly coincident 
with an optical/UltraVista/IRAC/radio source (separa- 
tion=0.09"), however it is within the sidelobe region of the 
brighter C0SLA-17-S source. Hence further follow-up with 
more complete uv-coverage is required to affirm this source. 
We find a photometric redshift of Zphot = 3.37 
sistent with the mm-to-radio based redshift of z, 
3.27 



-0.60 
-0.49- 



-1-0.14 
-0.22' 

mm/radio 



This SMG is detected by both LABOCA and 
AzTEC/ASTE surveys, and the flux ratio using these two 
surveys suggests a spectral index of 2.08. Using this value to 
scale the combined 1.3 mm PdBI fluxes to 870 fim (1.1 mm) 
we find a fiux density of 15.2 + 2.2 mJy (9.2 + 1.4 mJy), 
consistent with the deboosted LABOCA (AzTEC/ASTE) 
fiux of 12.5 + 3.2 mJy (7.51};? mJy). 

COSLA-18 

COSLA-18 is detected at S/N = 4.5 {a = 10 00 43.19, 
6 = +02 05 19.17). Its 1.3 mm fiux density is Fi.3mm = 
2.15 ± 0.48 mJy. Scaled to 870 /xm (1.1 mm) this yields 
a fiux density of 7.2 ± 1.6 mJy (3.5 + 0.8 mJy), consis- 
tent with the deboosted LABOCA (AzTEC/ASTE) fiuxes 
of 10.0 ± 2.6 mJy {3.8ti^ mJy). The source is coincident 
with an optical/UltraVista/IRAC/radio source (separation 
= 0.67"; 1+ = 28.96). Our photometric redshift compu- 
tation yields 2 almost equally probable photometric red- 
shifts at z ^ 3 and z 5. The significant radio detection 
(^'i.4GHz = 60.1 ± 8.9 fiJy) would argue in favor of the 
lower redshift solution, consistent with the mm-to-radio 
flux based redshift of ^mm/radio — 2.40lg;34. Thus, here 
we adopt the low-redshift solution for this source, yield- 
ing Zphot = 2.90^0 43, noting that a second solution of 

^phot 



4.92+°-^^ is possible. 



-0.34 
COSLA-19 

Cosla-19 is detected at S/N = 4.1 with a 1.3 mm flux 
density of 3.17 ± 0.76 mJy. Scaling this flux to 1.2 mm, 
1.1 mm, and 870 /iJy yields fluxes of 4.1 ± 1.0 mJy, 
5.1 ± 1.2 mJy, 10.7+ 2.6 mJy, respectively, consistent with 
the deduced MAMBO, AzTEC, and LABOCA fluxes of 
5.55 + 0.9 mJy, 5.3^1^2 niJy, and 7.4+1.8 mJy, respectively. 
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The closest multi-wavelength source to Cosla-19 is an opti- 
cal/UltraVista source 2.0" away (i+ = 25.66). Such a sepa- 
ration makes it unlikely that this source is the counterpart 
of the mm-detection although given the mm-resolution and 
significance a ~ 0.8" positional uncertainty is expected. A 
2(7 = 16.1 /iJy radio peak is associated with the PdBI mm 
peak yielding a mm-to-radio-flux ratio based redshift for 

Cosla-19 of Znun/radio 



optical/UltraVista/IRAC/radio source (separation= 0.55"; 

— 27.24). We find a photometric redshift of 2:pho' ~ 
1.91 



+ 1.75 
^0.64- 



o qS;+l-62 
•^■^°-0.90- 



COSLA-23 

Within the COSLA-23 LABOCA beam two significant 
(S/N > 5) sources are found in the PdBI 1.3 mm map. 
COSLA-23-N is detected at S/N = 7.3 at a 10 00 
10.161, S = +02 13 34.95. It is coincident with an op- 
tical/UltraVista/IRAC/radio source (separation^ 0.44"; 
i+ = 26.3). COSLA-23-S is detected at S/N = 6.2 at a = 10 
00 10.070, S = +02 13 26.87. It can be matched to an op- 
tical/IRAC source (separation= 0.87", i+ = 28.49), but 
it is not detected in the radio map. Based on the multi- 
wavelength photometry of the counterparts we find photo- 
metric redshifts of Zphot = 4.00^;^;^^ (COSLA-23-N) and 
2phot = ^.bStlil (COSLA-23-S). 

The 1.3 mm flux densities for COSLA-23-N and 
COSLA-23-S are 3.42 ± 0.47 mJy, and 3.70 ± 0.60 mJy, 
respectively. Only COSLA-23-N is within the MAMBO 11" 
beam, and the scaled 1.3 mm flux (4.4± 0.6 mJy) agrees 
weU with the COSBO-2 flux (5.77±0.9 mJy). 

COSLA-25 

No significant source is present in the 1.3 mm map. No 
statistical counterpart could be associated with this SMG. 
Therefore, it remains unclear whether the LABOCA source 
is spurious, breaks up into multiple components at 1.5" 
angular resolution, or simply is below the PdBI detection 
limit. 

COSLA-30 

No significant source is present in the 1.3 mm map. No 
statistical counterpart could be associated with this SMG. 
Therefore, it remains unclear whether the LABOCA source 
is spurious, breaks up into multiple components at 1.5" 
angular resolution, or simply is below the PdBI detection 
limit. 

COSLA-33 

The most prominent feature in the PdBI map within 
the LABOCA beam is a 3.1cr peak 6.0" away from 
the LABOCA source center that can be associated 
with an optical/UltraVista/IRAC source (i+ ~ 25.2, 
separation=0.95"). Its 1.3 mm flux density is 1.78 ± 
0.58 mJy which scales to 6.02 ± 1.95 mJy at 870 /xm, 
in good agreement with the deboosted LABOCA flux 
(6.8 ± 1.1 mJy). Given the low significance of the 1.3 mm 
source further follow-up is required to affirm its reality. 

COSLA-35 

COSLA-35 is detected at a signal-to-noise of S/N = 

4.2 (a = 10 00 23.65, S = +02 21 55.22). Its 1.3 mm 
flux density is Fi.3i„m — 2.15 ±0.51 mJy. This flux 
scaled to 870 fim (1.1 mm) yields a flux density of 

7.3 ±1.7 mJy (3.5 ± 0.8 mJy), consistent with the ob- 
served LABOCA (AzTEC/ASTE) fluxes of 8.2 ± 1.1 mJy 
(S.l^x'^ mJy). The mm-detection is coincident with an 



COSLA-38 

COSLA-38 is detected at S/N = 4.4 at a = 10 00 
12.59, S = +02 14 44.31, 14.8" away from the LABOCA 
source position (thus essentially outside the LABOCA 
beam; FWHM=27"). It is however only 0.67" away from 
the MAMBO source Cosbo-19, and coincident with a 
radio/UltraVista/IRAC/optical source (separation=0.23" , 
i+ — 24.08). We infer a 1.3 mm flux density of Fi.smm — 
8.19± 1.85 mJy, which should however be treated with cau- 
tion as the correction for the primary beam response at 
that distance from the PdBI phase center applied to the 
flux is about a factor of 4. We find a photometric redshift 
of zphot = 2.U-toll for this SMG. 

COSLA-40 

COSLA-40 is detected at S/N = 3.4 (a = 09 59 
25.91, S = +02 19 56.40), 11.3" away from the LABOCA 
source center. Its 1.3 mm fiux density is Fi.smm = 3.41 ± 

1.02 mJy. Scaling this flux to 870 fim yields a flux den- 
sity of 11.5 ± 3.4 mJy, in very good agreement with the 
deboosted LABOCA flux (Fsrof^m = H-l ± 3.4 mJy). 
The source is coincident with an optical source (separa- 
tion=0.51"; i+ = 25.52), but not detected in the radio. 
Given the expected flux density, and the coincidence of the 
source with an optical detection we assume this source to 
be real. We find a photometric redshift of 2;phot — l-30to^ii, 
but we note that dips are also found at lower and higher 
redshift values, and that the mm-to-radio flux based red- 
shift suggests Zmm/radio ^ 4.5. 

COSLA-47 

COSLA-47 is detected at a = 10 00 33.350, S = +02 26 
01.66 and S/N = 5.3, 6.4" away from the LABOCA source 
center. Its 1.3 mm flux density is 3.11 ± 0.59 mJy, and con- 
sistent with the LABOCA/AzTEC fluxes, when scaled to 
these frequencies. The PdBI source is coincident (separa- 
tion=0.48") with a source independently detected at opti- 
cal, IR, and radio wavelengths. We find a well constrained 
photometric redshift of Zphot = 2.36jlQ'24. 

Within the LABOCA beam several more S/N > 4 peaks 
can be associated with optical/UltraVista/IRAC sources. 
They are however within sidelobe contaminated regions. 
This LABOCA SMG may be a blend of several sources, 
but further follow-up is required to confirm this. 

COSLA-48 

No significant source is present in the 1. 3 mm map. The 
statistical counterpart association (see Sec. 5.2.1 and Fig. [5] 
for details) suggests two potential tentative counterparts 
to this LABOCA source. The sum of the extracted 1.3 mm 
fluxes (taken as maximum flux within a circular area of 
1" in radius centered at the statistical counterpart), cor- 
rected for the primary beam response, is 1.56 mJy. This 
flux, scaled to 870 ^m (5.2 mJy) is in very good agree- 
ment with the LABOCA flux of 6.1 ± 1.7 mJy (see Fig. [2]). 
This suggests that the LABOCA source may be fainter at 

1.3 mm than can be detected given our PdBI sensitivity 
and it breaks up into multiple components when observed 
at 1.5" resolution. 
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COSLA-50 

No significant source is present in the 1. 3 mm map. The 
statistical counterpart association (see Sec. |5. 2.1 and Fig.js] 
for details) suggests two potential (robust and tentative) 
counterparts to this LADOGA source. The sum of the ex- 
tracted 1.3 mm fluxes (taken as maximum flux within a cir- 
cular area of 1" in radius centered at the statistical counter- 
part), corrected for the primary beam response, is 2.61 mJy. 
When scaled to 870 fim this flux (8.7 mJy) is fairly consis- 
tent with the LABOCA flux of 5.6 ± 1.6 mJy (see Fig. ^. 
This suggests that the LABOCA source may be fainter at 
1.3 mm than can be detected given our PdBI sensitivity 
and it breaks up into multiple components when observed 
at 1.5" resolution. 

COSLA-51 

No significant source is present in the 1. 3 mm map. The 
statistical counterpart association (see Sec. 5.2.1 and Fig. [5] 
for details) suggests one robust potential counterpart to 
this LABOCA source. The extracted 1.3 mm flux (taken as 
maximum flux within a circular area of 1" in radius centered 
at the statistical counterpart), corrected for the primary 
beam response, is 1.27 mJy, consistent with the LABOCA 
flux of 6.2 ± 1.7 mJy when scaled to 870 /im (4.5 mJy; see 
Fig. [2]). This suggests that the LABOCA source is fainter 
at 1.3 mm than can be detected given our PdBI sensitivity 
and/or it breaks up into multiple components when ob- 
served at 1.5" resolution. 

COSLA-54 

COSLA-54 is detected at S/N = 5.0 {a =09 58 37.99, 
S = +02 14 08.52), 7.6" away from the LABOCA source 
center. Its 1.3 mm flux density is i^i.smm = 3.26±0.65 mJy. 
Scaling this flux to 870 /im (1.1 mm) yields a flux density 
of 11.0 ± 2.2 mJy (5.3 ± 1.1 mJy), in agreement with the 
deboosted LABOCA (AzTEC/ASTE) flux of 11.6±4.1 mJy 
mJy). The mm-detection can be associated with an 
optical/IRAC/radio source (separation= 0.75"; i+=25.21). 
We find a photometric redshift of Zphot = 2.64;to 26- 

COSLA-62 

No significant source is present in the 1.3 mm map. No 
statistical counterpart could be associated with this SMG. 
Therefore, it remains unclear whether the LABOCA source 
is spurious, breaks up into multiple components at 1.5" 
angular resolution, or simply is below the PdBI detection 
limit. 

COSLA-128 

COSLA-128 is detected at S/N = 4.8 (a = 10 01 37.99, 
S = +02 23 26.50). Its 1.3 mm flux density is Fi.gmm = 
4.50±0.94 mJy. Scaling this flux to 870 fim (1.1 mm) yields 
a flux density of 15.2±3.2 mJy (7.3±1.5 mJy), in agreement 
with the LABOCA (AzTEC/ASTE) flux of 11.0 ±3.5 mJy 
(4.4 ± 1.1 mJy). The source is coincident with an optical 
detection (no MIR/radio; separation= 0.55"; i+ = 26.57). 
We find a photometric redshift of Zphot = O.IO''""'^^ 
secondary and tertiary possible solutions at z • 

2-3. 



0.00' with 
' 1.2, and 



COSLA-161 

COSLA-161 is detected at S/N = 3.5 (a = 10 00 16.150, 
5 = +02 12 38.27). Its 1.3 mm flux density is Fi.3„u„ = 



2.54 ± 0.674 mJy. Scaling this flux to 870 /im, 1.1 mm, 
and 1.2 mm, using a spectral index of 3, yields 10.1 ± 4.8, 
4.1±1.2, and 3.3±0.9 mJy, respectively. This is in very good 
agreement with the deboosted LABOCA (10.1 ± 4.8mJy), 
AzTEC/ASTE (3.2 + 1.1 mJy), and MAMBO fluxes (1.4 + 
0.9 mJy). 

The SMG is coincident with an optical/IR/radio source 
with an available (VIMOS/IMACS) spectrum at ^spec = 
0.187. The source is also detected by Chandra in the X-rays, 
and we find a 0.5-2 keV band flux of 1.9 + 0.8 erg s~^ cm~^. 
At a redshift of 0.187 this corresponds to a bolometric X- 
ray luminosity (0.1-10 keV) of (6.2 + 2.8) x 10^° erg s'^ 
(assuming a power law spectral shape with photon index 
1.8). Given this X-ray luminosity it is not clear whether it 
arises from star-formation processes or a low-power AGN. 

It is interesting that a second radio source is present 
within the LABOCA beam (zgpec = 2.947), and is not as- 
sociated with mm-emission (however there is a 2.3ct peak 
at its position in the PdBI map). 



Appendix B: Notes on the 1.1 mm-selected sample 

Our l.lmm-selected sample is based on the SMA follow- 
up of 15 brightest SMGs drawn from the 1.1 mm 
AzTEC/JCMT-COSMOS survey at 18" angular resolu- 
tion (AzTEC-1 to AzTEC-15; see tabiinterf ; Younger 
et al. 2007, 2009). Detailed notes on individual targets 
are given in Younger et al. (2007, 2009). Here we have 
extracted the multi-wavelength photometry, tabulated in 
tab:aztecphot , for the counterparts of these SMGs us- 
ing the deep COSMOS multi-wavelength catalog, with 
UltraVista data added. The photometry in the I RAC bands 
had to be deblended for AzTEC-8 (see Fig. 



B.l| ), and 

that for AzTEC-10 had to specifically be extracted as this 
source was not present in the catalog (see Younger et al. 
2009). The photometry extraction and deblending were 
performed following the procedure described in detail by 
Smolcic et al. (2012). Furthermore, AzTEC-11 is a pecu- 
liar source that required particular attention. Younger et 
al. (2009) find that the SMA detection is best fit by a 
double Gaussian, suggesting a multiple component (N & 
S) source, labeled AzTEC-ll-N and AzTEC-ll-Sp] They 
present three positions for this SMG: i) AzTEC-11 when 
the SMA detection is fit using a single-Gaussian, and ii) 
AzTEC-ll-N and AzTEC-ll-S when the SMA detection is 
fit using a double Gaussian. AzTEC-11 is coincident with 
an optical/MIR/radio source with a spectroscopic redshift 
(2;spg(. = 1.599). AzTEC-ll-S (which is actually the north- 
ern component of the source) cannot be matched to a multi- 
wavelength counterpart in the deep COSMOS maps. Thus, 
given the rms in the 20 cm VLA-COSMOS survey we es- 



timate a mm-to-radio based redshift of 



[i/radi( 



> 2.58. 



AzTEC-ll-N (which is actually the southern component of 
the source) has an independent UltraVista and IRAC coun- 
terpart. To extract its photometry we have deblended the 
counterpart of AzTEC-ll-N by subtracting a 2D-Gaussian 
from the maps at the position of the counterpart of AzTEC- 
11, In Fig. |B.2| we show the deblended maps for AzTEC- 
11-N. We find a photometric redshift of ^phot = 1-51 
for this component. 



+0.41 
-0.92 



Note that the N 
Younger et al. 2009). 



S labels are inverted (see Tab. 1 in 
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Table A.l. Photometry table for our LABOCA SMGs with PdBI detections (magnitudes hsted are total AB magnitudes 
corrected for galactic extinction) 
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Table B.l. Photometry table for our l.lmm-selected sample (magnitudes listed are total AB magnitudes corrected for 
galactic extinction) 
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* AzTEC-2, at a spectroscopic redshift of 1.125 (Balokovic et al., in prep), is heavily blended 
by a bright, extended foreground galaxy. 
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Fig. B.l. Deblending of AzTEC-8 in Spitzer/IRAC bands. 



34 



V. Smolcic et al.: mm-imaging of COSMOS SMGs: Redshift distribution 



blended 



debiend«d 



J-b4nd; blended 



deblertded 




-s 



■ 















'■r . ■ ' ' ^ 

















5 -S 



S Q -S 



H-band: blended 



deblend«d 

T? — I — r 



Ks-bondr blended 



d«bler>ded 




S a -5 

*rt S»eondl 




5 0-5 
*jie S*cen4t 



! -5 



J.e^jfn; blended 



deb tended 




4.5|im; blended 
— ' I ■ ' I ■ I 



5.8fiiTi; blended 
■I '."J 



4 I 0-2-4 

debiended 





Fig. B.2. Deblending of AzTEC-ll-N in UltraVista YJHK and Spitzer/IRAC bands. 



